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Qh! Abstract 



The AcerMC Monte Carlo Event Generator is dedicated for the generation of 
Standard Model background processes at pp LHC collisions. The program itself 
provides a library of the massive matrix elements and phase space modules for 
generation of a set of selected processes: gg, qq — » tibb, qqW{-^ lv)bb, qqW{-^> lv)ti, 
99>QQ ~^ U)bb, gg,qq — > Zj^f{-^ tt,vv,W)tt and complete electroweak 

gg — ► (Z/W/'j* —>)bbti process. The hard process event, generated with one of these 
modules, can be completed by the initial and final state radiation, hadronisation 
and decays, simulated with either PYTHIA or HERWIG Monte Carlo event generator. 
Interfaces to both of these generators are provided in the distribution version. The 
matrix element codes have been derived with the help of the MADGRAPH package. 
The phase-space generation is based on the multi-channel self-optimising approach 
as proposed in NEXTCALIBUR event generator. Eventually, additional smoothing of 
the phase space was obtained by using a modified ac -VEGAS routine in order to 
improve the generation efficiency. 
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1 PROGRAM SUMMARY 



Title of the program: AcerMC version 1.0 

Operating system: Linux 

Programming language: FORTRAN 77 with popular extensions. 
External libraries: CERNLIB. 

Size of the compressed distribution directory: about 3.6 MB. The distribution 
includes modified versions of PYTHIA 6 . 2, HERWIG 6 . 3 and HELAS libraries. 
Key words: Standard Model backgrounds at LHC, massive matrix elements, 
Monte Carlo generator, heavy flavor production, multi-channel phase-space 
generation. 

Nature of physical problem: Despite a large repertoire of processes implemented 
for generation in event generators like PYTHIA [1] or HERWIG [2] a number of 
background processes, crucial for studying expected physics potential of the 
LHC experiments is missing. For some of these processes the matrix element 
expressions are rather lengthly and/or to achieve a reasonable generation ef- 
ficiency it is necessary to tailor the phase-space selection procedure to the 
dynamics of the process. That is why it is not practical to imagine that any of 
the above general purpose generators will contain every, or even only observ- 
able, processes which will occur at LHC collisions. A more practical solution 
can be found in a library of dedicated matrix-element-based generators, with 
the standardised interfaces like that proposed in [3] , to the more universal one 
which is used to complete the event generation. 

Method of solution: The AcerMC Event Generator provides itself library of 
the matrix-element-based generators for a few example processes. The initial- 
and final- state showers, beam remnants and underlying events, fragmenta- 
tion and remaining decays are supposed to be performed by the other uni- 
versal generator to which this one is interfaced. We will call it supervising 
generator. The interfaces to PYTHIA 6 . 2 and HERWIG 6.3, as such genera- 
tors, are provided. At present, the following matrix-element-based processes 
have been implemented: gg, qq — > tibb, qq — > W (— > tv)bb; qq — > W{— > £v)tt, 
gg,qq -> Z/Y(^ #)&&; gg,qq -> Z/Y{^ tt,vv,bb)tt and complete EW 
gg — > (Z/W/j* ->)tibb. Both interfaces allow the use of the PDFLIB library of 
parton density functions. 

Restriction on the complexity of the problem: The package is dedicated for the 
14 TeV pp collision simulated in the LHC environment. In particular, using it 
for the 2 TeV pp collision although technically possible might not be efficient 
and would require dedicated optimisation of the phase space generation. The 
consistency between results of the complete generation using PYTHIA 6 . 2 or 
HERWIG 6 . 3 interfaces is technically limited by the different approaches taken 
in both these generators for evaluating o:qcd and qlqed couplings and by the 
different models for fragmentation/hadronisation. For the consistency check, 
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in the AcerMC library contains native coded definitions of the olqcd and 
ocqed- Using these native definitions leads to the same total cross-sections 
both with PYTHIA 6 . 2 or HERWIG 6 . 3 interfaces. Typical running time: On 
an PHI 800 MHz PC it amounts to ~ 0.05 — > 1.1 events/sec, depending on 
the choice of process. 

[1]. T. Sjostrand et al., High energy physics generation with PYTHIA 6.2, 
eprint |hep-ph/0l08^ , LU-TP 01-21, August 2001. 



[2]. G. Marchesini et al, Comp. Phys. Commun. 67 (1992) 465, G. Corcella 
et al, JHEP 0101 (2001) 010. 



[3]. E. Boos at al., Generic user process interface for event generators, |hcp- 
ph/0109068 . 
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2 Introduction 



Despite a large repertoire of processes implemented for generation in the uni- 
versal generators like PYTHIA [1] or HERWIG [2] a number of Standard Model 
background processes for studying expected physics potential of the LHC ex- 
periments is still missing. For some of these processes the matrix element 
expressions are rather lengthy and/or to achieve a reasonable generation ef- 
ficiency it is necessary to tailor the phase-space selection procedure to the 
dynamics of the process. Due to this fact it cannot be expected that any of 
the universal purpose generators will contain every proces that is expected to 
participate at LHC pp collisions with an observable rate. A more practical so- 
lution could come in form of dedicated matrix-element-based generators with 
standardised interfaces, like the one proposed in [3], to the more general ones 
which are used to complete event generation. 

The AcerMC Monte Carlo Event Generator follows up on this idea. It is ded- 
icated for the simulation of the specific Standard Model background processes 
at LHC collisions: the gg, qq — > tibb, qq — > W(-^ £u)bb; qq — > W (— > £u)tt : 
gg,qq -> Z/Y(^ M)bb; gg,qq -> Z/Y(^ tl,vv,Vb)tt and complete EW 
gg — > (Z/W/Y —>)ttbb. They are characterised by the presence of the heavy 
flavour jets and multiple isolated leptons in the final state. For the Higgs boson 
searches, the ttH, ZH, WH with H -> 66, the gg -> H with H -> ZZ* -> U, 
the bbh/H/A with h/H/A — > rr, fj,fj, are the most obvious examples of signals 
where the implemented processes would contribute to the dominant irreducible 
backgrounds. The same background processes should also be considered for 
e.g. estimating the observability of SUSY events with a signature of multi-b-jet 
and multi-lepton production. 

The program itself provides library of the massive matrix elements and phase 
space modules for the generation of a few selected 2^4 processes. The 
hard process event, generated with these modules, can be completed by the 
initial and final state radiation, hadronisation and decays, simulated with 
either PYTHIA 6 . 2 or HERWIG 6 . 3 Monte Carlo Event Generators. These will 
subsequently be called the Supervising Generators. Interfaces of AcerMC to 
both, PYTHIA 6 . 2 or HERWIG 6 . 3 generators, are provided in the distribution 
version. The AcerMC also uses several other external libraries: CERNLIB, 
HELAS [4], VEGAS [5]. The matrix element codes have been derived with the help 
of MADGRAPH [6] package. The achieved typical efficiency for the generation of 
unweighted events is of 20% - 30%, rather high given a complicated topology 
of the implemented processes. 

The very first version of this library, interfaced to PYTHIA 6 . 1 within the 
standard of the so-called external processes (i.e. stand-alone implementations 
of hard processes interfaced to PYTHIA for further treatment of ISR/FSR and 
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hadronisation, c.f [1] ), was already available and used by ATLAS Collabo- 
ration for physics simulation studies since several months. The documenta- 
tion, for the processes implemented in this early version: qq — > W(— * £u)bb, 
gg,qq — > Z/Y(^ t£)bb and gg,qq — > tibb, can be found respectively in [7], 
[8], [9]. Since then, when upgrading to the AcerMC standard, the efficiency 
has been significantly improved thanks to the additional optimisation step in 
the phase space generation. Also, the interface standard was changed from 
PYTHIA 6.1 to PYTHIA 6 . 2 conventions, an interface to HERWIG 6 . 3 gen- 
erator was introduced and the native AcerMC calculations of the cxqed and 
olqcd couplings were coded to allow for consistent benchmarking between 
results obtained with PYTHIA and HERWIG as supervising generators. As a sig- 
nificant extension, the: qq — > W(— > £is)tt; gg,qq — > Zj^j*{— > ££,vu,bb)ti and 
complete electroweak gg — > (Z/W/Y —>)ttbb processes were added, which 
have been implemented for the first time in the AcerMC library. 

The outline of this paper is as follows. In Section 3, we describe physics mo- 
tivation for implementing each of the above processes and we collect some 
numerical results (plots, tables) which can be used as benchmarks. In Section 
4 we describe the overall Monte Carlo algorithm. Section 5 gives details on 
the structure of the program. Section 6 collects information on how to use this 
program and existing interfaces to PYTHIA 6 . 2 and HERWIG 6 . 3. Summary, 
Section 7, closes the paper. Appendix A documents sets of Feynman diagrams 
used for calculation of the matrix element for each subprocess, Appendices B 
and C give examples of the input/output of the program. 
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3 Physics content 



The physics programme of the general purpose LHC experiments, ATLAS [10] 
and CMS [11], focuses on the searches for the New Physics with the distinctive 
signatures indicating production of the Higgs boson, SUSY particles, exotic 
particles, etc. The expected environment will in most cases be very difficult, 
with the signal to background ratio being quite low, on the level of a few 
percent after final selection in the signal window. 

Efficient and reliable Monte Carlo generators, which allow one to understand 
and predict background contributions, are becoming the key point to the dis- 
covery. As the cross-section for signal events is rather low, even rare Standard 
Model processes might become the overwhelming background in such searches. 
In several cases, generation of such a process is not implemented in the general 
purpose Monte Carlo generators, when the complicated phase space behaviour 
requires dedicated (and often rather complex) pre-sampling, whilst the gen- 
eral purpose Monte Carlo generators due to a large number of implemented 
processes tend to use simpler (albeit more generic) phase space sampling algo- 
rithms. In addition, the matrix element for these processes is often lengthy and 
thus requiring complicated calculations. Only recently, with the appearance of 
modern techniques for automatic computations, their availability on demand 
became feasible for the tree-type processes. With the computation power be- 
coming more and more easily available, even very complicated formulas can 
now be calculated within a reasonable time frame. 

The physics processes implemented in AcerMC library represent such a set 
of cases. They are all being key background processes for the discovery in 
the channels characterised by the presence of the heavy flavour jets and/or 
multiple isolated leptons. For the Higgs boson searches, the tiH, ZH, WH with 
H -> bb, the gg -> H with H -> ZZ* -> U, the bbh/H/A with h/H/A -> 
tt, jiji are the most obvious examples of such channels. 

It is not always the case that the matrix element calculations in the lowest 
order for a given topology represent the total expected background of a given 
type. This particularly concerns the heavy flavour content of the event. The 
heavy flavour in a given event might occur in the hard process of a much 
simpler topology, as the effect of including higher order QCD corrections (eg. 
in the shower mechanism). This is the case for the b-quarks present in the 
inclusive Z-boson or W-boson production, which has a total cross-section or- 
ders of magnitude higher than the discussed matrix-element-based Wbb or Zbb 
production. Nevertheless, the matrix-element-based calculation is a very good 
reference point to compare with parton shower approaches in different frag- 
mentation/hadronisation models. It also helps to study matching procedures 
between calculations in a fixed oiqcd order and parton shower approaches. For 
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very exclusive hard topologies matrix-element-based calculations represent a 
much more conservative approximation than the parton shower ones [12]. 

Let us shortly discuss the motivation for these few Standard Model background 
processes which are implemented in the AcerMC 1.0 library. 

The tibb production at LHC is a dominant irreducible background for the 
Standard Model (SM) and Minimal Supersymmetric Standard Model (MSSM) 
Higgs boson search in the associated production, tiH, followed by the decay 
H — ► bb. The potential for the observability of this channel has been carefully 
studied and documented in [10] and [13]. Proposed analysis requires identify- 
ing four b-jets, reconstruction of both top-quarks in the hadronic and leptonic 
mode and visibility of the peak in the invariant mass distribution of the re- 
maining b-jets. The irreducible tibb background contributes about 60-70% of 
the total background from the ti events (tibb, tibj, tijj). 

The Wbb production at LHC is recognised as a substantial irreducible back- 
ground for the Standard Model (SM) and Minimal Supersymmetric Standard 
Model (MSSM) Higgs boson search in the associated production, WH, fol- 
lowed by the decay H — > bb. The massive matrix element for qq — > Wg*(—> bb) 
process has been calculated [14] and interfaced with HERWIG 5 . 6 Monte Carlo 
[15] already a few years ago. A more recent implementation of the Wbb + 
multi-jet final states is available from [16]. Recently, the massless matrix el- 
ement has been implemented in the general purpose Monte Carlo program 
MCFM [17], where the radiative corrections to this process are also addressed. 
Another implementation of the qq — > W(— > iv)g*(-^ bb) massive matrix el- 
ements, with the interface to PYTHIA 6.1 became available in [7]. The Ac- 
erMC library discussed here includes even more efficient implementation of 
the algorithm presented in [7]. 

The Wti production at LHC has to our knowledge not been implemented 
in the publicly available code so farQ. It is of interest because it contributes 
an overwhelming background [18] for the measurement of the Standard Model 
Higgs self-couplings at LHC in the most promising channel pp — > HH — > 
WWWW. 

The Z/~f*( — ► tt)bb production at LHC has since several years been recog- 
nised as one of the most substantial irreducible (or reducible) backgrounds 
for the several Standard Model (SM) and Minimal Supersymmetric Standard 
Model (MSSM) Higgs boson decay modes as well as for observability of the 
SUSY particles. There is a rather wide spectrum of regions of interest for this 
background. In all cases the leptonic Z/j* decay is asked for, but events with 
di-lepton invariant mass around the mass of the Z-boson mass or with the 

2 We thank M. L. Mangano for bringing this process to our attention and for 
providing benchmark numbers for verifying the total cross-section. 
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masses above or below the resonance peak could be of interest. The presented 
process enters an analysis either by the accompanying b-quarks being tagged 
as b-jets, or by the presence of leptons from the b-quark semi-leptonic decays 
in these events, in both cases thus contributing to the respective backgrounds. 

Good understanding of this background, and having a credible Monte Carlo 
generator available, which allows studying of expected acceptances for differ- 
ent final states topologies, is crucial. Despite a very large effort taken in time 
of the preparation of the Aachen Workshop [19], such well established Monte 
Carlo generator was missing for several years. The matrix element for the 
gg — > Zbb — > bbtt production has been published already in [20] and in time 
of Aachen Workshop implemented into EURO JET Monte Carlo [21]. But that 
generator was not giving the possibility for having fully generated hadronic 
event with modelled initial and final state radiation and hadronisation, for 
analyses presented in [22] it was interfaced to PYTHIA 5 . 6 [23] Monte Carlo. 
This is no longer supported and available at present. The same matrix ele- 
ment has been directly implemented into PYTHIA 5.7 [23]. However, with this 
implementation the algorithm for the phase space generation never working 
credibly, it has finally been removed from the version PYTHIA 6 . 1 [24]. A year 
ago, the massless matrix elements for gg,qq — > Zbb processes have been im- 
plemented in the general purpose Monte Carlo program MCFM [17]. In that 
implementation radiative corrections to this process are addressed as well. 
The massive matrix element implementation is also present in the very recent 
version of HERWIG 6.3 [2]. At that time the gg,qq — > Zj^fbb — > ttbb massive 
matrix elements, with the interface to PYTHIA 6.1 became available [8]. The 
AcerMC library discussed here includes more efficient implementation of the 
algorithm presented in [8]. 

The Z/*y*( — > ££,isv,bb)ti production at LHC is an irreducible back- 
ground to the Higgs search in the invisible decay mode (case of Z — > vv) in the 
production with association to the top-quark pair [25]. With the Zj^f{— > bb) 
it is also an irreducible resonant background to the Higgs search in the tiH 
production channel but with the Higgs boson decaying to the b-quark pair 
[13]. 

The complete EW production of the gg — > (Z/W/~y* — final state 
is also provided. It can be considered as a benchmark for the previous process, 
where only the diagrams with resonant gg — > (Z/^y* —>)bbti are included. It 
thus allows the verification of the question, whether the EW resonant con- 
tribution is sufficient in case of studying the ttbb background away from the 
Z-boson peak, like for the tiH with Higgs-boson mass of 120 GeV. 

This completes the list of the native AcerMC processes implemented so far. 
Having all these different production processes implemented in the consistent 
framework, which can also be directly used for generating standard subpro- 
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cesses implemented in either PYTHIA or HERWIG Monte Carlo, represents a very 
convenient environment xor several phenomenological studies dedicated to the 
LHC physics. 

For the cases, where radiative photon emission from final state leptons is im- 
portant the package PHOTOS [26] can be used in the chain of event generation. 
In similar way also package TAUOLA [27] can be interfaced directly to the gen- 
eration chain and used for events generation in cases where more detailed 
treatment of the tau-lepton decay and including spin correlations effects is 
relevant. 

In the following subsections we discuss in more detail implementation of each 
subprocess. We also give benchmark Tables with the total cross-sections ob- 
tained with AcerMC processes but different implementations and setting of 
«qcd {Qqcd)'- the native AcerMC, PYTHIA and HERWIG ones. For a more de- 
tailed discussion on this topic the reader is referred to Section 4.3. If the 
native AcerMC definition is used, the same cross-section is obtained either 
with PYTHIA or HERWIG generation chains. 

3.1 The gg, qq — > tibb processes 

In the implementation discussed here, the matrix element was derived using 
the MADGRAPH package [6]. These matrix elements are not covering the decay of 
the top-quarks, the latter are considered as massive final states of the process. 
The top-quark decays is than performed by the supervising generator. Rather 
important spin effects (spin correlations) in the top decays are therefore not 
yet included. The similar solution, like for tau decay in the Z-boson production 
process discussed in [28], is planned to be implemented here in the near future. 

As a benchmark, the processes gg, qq — > tibb have been simulated for pp col- 
lisions with 14 TeV centre-of-mass energy and CTEQ5L [29] parton density 
functions, using event generation with massive 2^4 matrix element imple- 
mented as an external process to PYTHIA 6 . 2 (see Section 4 and 5). The decays 
of the top-quarks have been left under control of PYTHIA 6 . 2 generator. The 
qq — > tibb subprocess contributes less than 10% of the total cross-section. 

The total cross-section is very sensitive to the choice of the QCD energy scale 
used for calculation of that process, thus indicating potentially large contribu- 
tions from higher order corrections. The same definition for the factorisation 
and renormalisation scale is used. The example values of the total cross-section 
for implemented choices of the QCD energy scale are given in Table 1. 

As a cross-check, the processes gg, qq — > tibb have been coded independently 
using the COMPHEP package [30]. The same set of diagrams was selected and 
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Table 1 

AcerMC cross-sections for the gg, qq — ► ttbb production at different choices of the 
QCD energy scale and uqcd implementations. The 14 TeV centre- of-mass energy 
and CTEQ5L parton density functions were used for the simulation with interfaces 
to PYTHIA 6.2 and HERWIG 6.3. The m H = 120 GeV and m t = 175 GeV were 
used for calculating the Qqcd ^ n ^ as ^ row °f ^is table. The default settings of 
a QCD as implemented in AcerMC, PYTHIA 6.2 and HERWIG 6.3 were used. 



Factorisation scale 


oqcd (1L) 
native AcerMC 


oqcd (1L) 
as in PYTHIA 6 . 2 


a QCD (2L) 
as in HERWIG 6 . 3 




gg — > ttbb 


Qqcd = S 


4.2 [pb] 


3.9 [pb] 


2.7 [pb] 


Qqcd = E(Pt 2 + ™?)/ 4 


10.3 [pb] 


10.2 [pb] 


6.4 [pb] 


Qqcd = E(Pt 2 )/ 4 


17.0 [pb] 


16.9 [pb] 


10.1 [pb] 


Qqcd = ( mt + m »/ 2 ) 2 


8.2 [pb] 


8.1 [pb] 


5.2 [pb] 




qq — > ttbb 


Qqcd = * 


0.30 [pb] 


0.29 [pb] 


0.20 [pb] 


Qqcd = I>t 2 + ™?)/ 4 


0.61 [pb] 


0.60 [pb] 


0.38 [pb] 


Qqcd = E^')/ 4 


0.91 [pb] 


0.90 [pb] 


0.54 [pb] 


Qqcd = ( mt + m »/ 2 ) 2 


0.52 [pb] 


0.51 [pb] 


0.33 [pb] 



only the integrating part of the package was used to calculate total cross- 
section. The choices for the QCD energy scale were kept consistent. A very 
good agreement between the cross-sections obtained with two independent 
calculation streams prepared for this study has been achieved [9]. 

One can observe a very strong scale dependence of the cross-section for the 
gg,qq — > ttbb process (c.f. Table 1). Factor four (!!) can be expected on the 
predicted cross-section when changing from the scale Qqcd = s to the scale 
Qqcd = < Pt >• This very strong dependence on the energy scale is 
also observed in the case of the tiH production, for recent discussion see [31]. 
There, the recommended central factorisation and renormalisation energy scale 
is fiQ = (mt + m#/2). Having in mind that the primary interest of evaluating 
this background is the Higgs search in the tiH production, i.e. with the fa- 
quark system being produced with the invariant mass of the expected Higgs 
boson, we have also introduced this central energy scale, with m H = 120 GeV 
as one of the possible choices. 

Fig. 1 shows the distributions of the Qqcd = \JQqcd (distributions have 
been normalised to one) for the ttbb events with the invariant mass of the b- 
jets system, calculated using the default PYTHIA (LO) oqcd implement at ionQ. 
rribb-jets = 120 ± 30 GeV. The distribution is well collimated around the 



3 This would makes distributions directly relevant for the ttH analysis. For details 
on the jet reconstruction see [32]. 



13 



average value when Qqcd * s defined as < > or < p\ > while it is much 
broader when Qqcd * s defined as s. The kinematic distributions are very sim- 
ilar in shape for separate gg — > tibb and qq — > tibb contributions. The total 
cross-section for accepted events as a function of the averaged Qq C d (f° r these 
events) is shown in the bottom plot. It can be noted that the cross-section 
decreases rather fast with the increasing value of the average < Qq C d >• 
Also shown is the ci a s {Qqcd) dependence scaled to match the cross-section at 
Qqcd = + m t ) with mn = 120 GeV, it being the only calculated 

cross-section point with a fixed scale. The other cross-sections are shown to 
follow the expected ol^{Qqcd) dependence rather well, while the deviations 
are induced by the parton density function dependence on the Qqcd scale, 
most notably at Qqcd = & value. The deviations induced by the parton den- 
sity functions dependence on the Qqcd sca le are different for the gg and qq 
contributions, as can be concluded from results given in Table 1. 

The series of plots illustrating the most relevant differential distributions for 
the top-quarks and b-quarks can be found in [9]. 




Fig. 1. Top: the Qq C d distributions for tibb events with the invariant mass of the 
b-jets system rribb-jets = 120 ± 30 GeV. Bottom: the total cross-section of accepted 
events as a function of the averaged Qq C d (f or these events). 
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3.2 The qq — s- VF(— ► £u)g*(^ bb) process 

The matrix element for the implemented process was again coded by using 
the MADGRAPH package [6]. This process is represented by only two Feynman 
diagrams, with quark exchange in the t-channel, leading to the production 
of the VT-boson and virtual gluon splitting into bb pair. Only the u, d, s, c 
quarks were considered in this implementation, the possibility of the b-quark 
in the initial state was omitted as expected to be negligible numerically (e.g.. 
\Vbc/V u d\ 2 ~ 0.002) but leading to several additional diagrams which would 
have to be included. The massive matrix element takes into account spin 
correlations in the W-boson decay and angular correlations between leptons 
and quarks. Due to the massive treatment of the final state fermions the 
amplitude has no singularities; the total cross-section is well defined. The 
effect from the H^-boson natural width and the P^-boson propagator are also 
properly included. 

As a benchmark, the process qq — > W(— » £v)g*(—> bb) has been simulated 
for pp collision with 14 TeV centre-of-mass energy. The total cross-section, 
including branching ratio for W — > Iv (single flavour) is 36.5 pb (CTEQ5L 
parton density functions, Q 2 = Myy, PYTHIA 6.2 interface)^ 

Table 2 

AcerMC production cross-sections for the qq — * Wbb with W — * ev decay (single 
flavour). The 14 TeV centre-of-mass energy and CTEQ5L parton density functions 
were used with different definitions of olqed, olqcd ( as in default PYTHIA 6.2 
and HERWIG 6.3) and several choices of the factorisation scale, olqed a nd olqcd 
implementations. 



Factorisation scale 


&QED, Q-QCD (1L) 
native AcerMC 


a QED, OLQCD (1L) 
as in PYTHIA 6 . 2 


OLQED, olqcd (2L) 
as in HERWIG 6 . 3 


Q Z = M 2 W 


36.5 [pb] 


36.4 [pb] 


29.5 [pb] 


Q 2 = 'a 


44.1 [pb] 


44.0 [pb] 


34.8 [pb] 




36.0 [pb] 


36.0 [pb] 


29.1 [pb] 


Q 2 = (*w + s l- b )/ 2 + P T w 


37.2 [pb] 


37.1 [pb] 


30.0 [pb] 



The dependence on the choice of the factorisation scale is rather modest (c.f. 
Table 2) and does not exceed 20% for the choices implemented in AcerMC 

4 This can be compared with the matrix element implementation to HERWIG 5.6, 
used in [33], [34], where originally this cross-section was estimated to 19.8 pb 
(CTEQ2L parton density functions) but, when implementing CTEQ5L parton den- 
sity functions and setting kinematic parameters to be in approximate accordance 
with PYTHIA defaults, rises to 36.0 pb, which is consistent with the AcerMC imple- 
mentation by taking into account the remaining differences in the two calculations 
(e.g. the former implementation uses an on-shell W boson in the ME calculation). 
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library. The variation of the cross-section due to different olqed and olqqd 
implementations and default settings is again evident; as one can expect the 
two- loop oiqcd implementation given in HERWIG gives a ~20 % lower cross- 
section when compared to the cases when native AcerMC and PYTHIA one- 
loop (xqcd were usedQ. 

The differential distributions of the qq — *> Wbb events turn out to be interest- 
ing when compared to the corresponding ones of the qq — *> Zbb and gg — *> Zbb 
events (generated with pure Z-boson exchange). Such comparison is well doc- 
umented in [8]. 

3.3 The qq ^W{-^ tv)g*{-^ tt) process 

The matrix elements, coded by the MADGRAPH package [6], are not covering the 
decay of the top-quarks; the latter are considered as massive final states of the 
process. The top decay is than performed by the supervising generator. As in 
the case of gg, qq — > ttbb process spin effects in the top decays are therefore not 
yet included. This process, although rare, contributes an overwhelming irre- 
ducible background to possible measurement of the Higgs-boson self-coupling 
in the EE -> WWWW decay mode [18]. 

Table 3 shows the expected AcerMC cross-sections for different choices of 
the energy scale and coupling [olqed, a QCo) definitions. One should notice 
the effect of almost a factor two from different choices of the energy scale. 

Table 3 

AcerMC 'production cross-sections for the qq — > Wti with primary W — > ev 
decay (single flavour). The 14 TeV centre- of -mass energy, CTEQ5L parton density 
functions with different factorisation scales and different definitions of the olqed 
and olqcd were used in the matrix element calculations. 



Factorisation scale 


OiQED, olqcd (1L) 
native AcerMC 


a QED, a QCD (1L) 
as in PYTHIA 6 . 2 


OLQED, OLQCD (2L) 
as in HERWIG 6 . 3 


fft — Af2 
^QCD ~ 1V1 W 


69.3 [fb] 


69.1 [fb] 


56.0 [fb] 


Qqcd = s ti 


40.9 [fb] 


39.9 [fb] 


33.9 [fb] 


Qqcd = M w+P T w 


59.7 [fb] 


59.5 [fb] 


48.8 [fb] 


Qqcd = ( s w + s 'u)/ 2 + P T w 


43.7[fb] 


42.8 [fb] 


36.0 [fb] 



5 While performing further comparisons of native AcerMC and PYTHIA processes 
we discovered a misinterpretation of our CKM matrix implementation. This correc- 
tion efectively changes the cross-section for qq — > Wbb and qq — * Wti processes by 
~10% compared to the draft versions of this paper, which is nevertheless still well 
within the physics precision of the program. The affected tables in this paper are 
already updated. 
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3.4 The gg, qq — > Z/j*(— > ii)bb processes 

The matrix elements, derived using the MADGRAPH package [6], properly take 
into account spin correlations in the Z-boson decay and angular correlations 
between leptons and quarks. Thank to keeping non-zero b-quark masses the 
amplitude has no singularities; the total cross-section is well defined. 

The full Z/Y exchange proves to be important: For events well below the 
Z-boson resonance the contribution from 7* becomes dominant; the 7* contri- 
bution is also sizeable in the high mass tail and increases proportionally with 
the effective mass of the di-lepton system. 

As a benchmark result, the process has been simulated for pp collisions at 
14 TeV centre-of-mass energy. The total cross-sections, including the branch- 
ing ratio for Z/'j* — » ii (single flavour) are given in Table 4 for different 
definitions of cyqed, olqcd couplings. 

Several differential benchmark distributions for leptons and b-quarks originat- 
ing from the hard process has been collected and discussed in [7]. 

Table 4 

AcerMC production cross-sections for the gg,qq — ► Zj^fbb with Z/7* — > ee 
decay (single flavour). The 14 TeV centre-of-mass energy, CTEQ5L parton den- 
sity functions and different definitions for the olqed, olqcd ( a $ in default PYTHIA 
6 . 2 and HERWIG 6.3) were used in the matrix element calculations. The threshold 
met > 10 GeV was used in the event generation. 



Factorisation scale 


otQED, oiqcd{1L) 
native AcerMC 


(XQED, CiQCD{lL) 

as in PYTHIA 6.2 


(XQED, CtQCD&L) 

as in HERWIG 6 . 3 




gg - Zh*bb 


Q 2 = m| 


26.4 [pb] 


26.4 [pb] 


20.5 [pb] 




qq — > Z/~f*bb 


Q 2 = ml 


4.3 [pb] 


4.3 [pb] 


3.3 [pb] 



3.5 The gg, qq — > Z/j*(— > ££, vv, bb)tt processes 

This process, in spite of having a very small cross-section at LHC energies, 
contributes as irreducible background to the ttH production at low masses. In 
case the Higgs boson is searched within the H — > bb mode, this contribution 
becomes less and less important with the Higgs boson mass moving away from 
the Z-boson mass. In case of the Higgs-boson search in the invisible decaying 
mode, the Z —> vv might be more relevant also for the higher masses, as the 
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mass peak cannot be reconstructed for signal events. The Z/'-f* — > ££ decay is 
of less interest, as the expected observability at LHC is very low (c.f. Table 
5). 

Table 5 

AcerMC production cross-sections for the gg,qq — > Z/^*ti with Z — > zv e zv e decay 
(single flavour). The 14 TeV centre- of -mass energy, CTEQ5L parton density func- 
tions and different definitions for the uqed, oiqcd (as in native AcerMC, default 
PYTHIA 6.2 and HERWIG 6.3) were used. The threshold m UeUe > 60 GeV was used 
in the event generation. 



Factorisation scale 


OiQED, CtQCD^L) 

native AcerMC 


CtQED, CtQCD^-L) 

as in PYTHIA 6 . 2 


CtQED, CXQCd{2L) 

as in HERWIG 6.3 




99 Z/j*(^ v e v e )tt 


Q-m z 


41.3 [pb] 


41.3 [pb] 


32.1 [pb] 




qq -► Z/j*{^ v e v e )tt 


Q 2 = m 2 z 


21.2 [pb] 


21.2 [pb] 


16.5 [pb] 



3.6 The electron) eak gg — > (Z/W/j* — >)i>b££ process 



One should be well aware, that the gg,qq — > Z/j*ti with Z/7* — > 66 does 
not represent a complete electroweak production of the tibb final state. Conse- 
quently, a separate implementation for generation of the complete set of such 
diagrams (including as well W-boson exchange) was addressed. In fact this 
final state leads to complicated pattern of the 72 Feynman diagrams. 

The contribution from all non-resonant channels is a dominant one for the 
inclusive cross-section, see Table 6. An almost factor 10 higher cross-section 
is calculated with the full electroweak gg — > {ZjWj^f —*)bbtt with respect to 
calculated with the gg — > {Z/'-f* — > bb)tt process only. The interesting region 
for the background estimates to the Higgs searches is the one with the mass of 
the bb system around 120 GeV, see Fig. 2. This contribution was not included 
in results discussed in [10], [13], so these analyses will require revisiting. One 
should also note that the electroweak gg — > (ZjWj^* —>)bbti inclusive cross- 
section is on the level of 10% of the QCD gg — > bbtt cross-section, see Table 1, 
for the same choice of the energy scale. But in the mass range around 120 GeV 
it is on the level of 50% of the QCD contribution, as clearly indicated in Fig. 2. 

The implementation of the EW qq — > (Z/W/ r y* —>)bbti is still lacking but will 
be added in the near future. One could expected, on the base of cross-sections 
in Table 5, that the contribution from the quark-antiquark anihilation will be 
also significant. 
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Table 6 

AcerMC froduction cross-sections for the electroweak gg — ► (Z/W/j* -^)bbti and 
only for the resonant gg — > (Z/j* — > T/ie TeF centre- of-mass energy 

and CTEQ5L parton density functions were used along with different definitions for 
the ctQED, olqcd (as in native AcerMC, default PYTHIA 6.2 and HERWIG 6.3). 
The threshold m b i > 10 Ge V was used for generation in the resonant case. The 
tuh = 120 Ge V was used for caluculation of the energy scale. 



Factorisation scale 


&QED, ciqcd(^L) 
native AcerMC 


O-QED, CHQCdO-L) 
as in PYTHIA 6 . 2 


CtQED, ocqcd(2L) 
as in HERWIG 6 . 3 




gg -> {Z/Wh* ->)bbtt 


Qqcd = * 


0.58 [pb] 


0.56 [pb] 


0.71 [pb] 


Qqcb = E(Pt 2 +™?)/4 


1.10 [pb] 


1.05 [pb] 


0.84 [pb] 


QqCD = E(fT 2 )/ 4 


1.50 [pb] 


1.50 [pb] 


1.16 [pb] 


Qqcd = (m t +m H /2) 2 


0.90 [pb] 


0.89 [pb] 


0.71 [pb] 




gg -» (Zh* -» 66)« 


O 2 — M 2 
^QCD ~ 1V1 Z 


0.092 [pb] 


0.092 [pb] 


0.071 [pb] 




Fig. 2. The invariant mass distribution of the b-quark system: Left: EW 
gg — ► (Z/W/j* —*)ttbb (solid) and only resonant gg — > (Z/j* — > 66) it (dashed); 
Right: 55 -> (Z/W/j* -^)tt66 (so/idj and QCD 55 -> (dashed). The en- 
ergy scale Qqcd = ( m * + m #/2) 2 with ran = 120 GeV was used for generating 
QCD and full electro-weak production. 



19 



4 Monte Carlo algorithm 

The conceptual motivation leading to the present implementation of AcerMC 
was to exploit the possibility of dedicated matrix-element-based generation in- 
terfaced to a more general event generator, called supervising event generator, 
which is subsequently used to complete the event generation procedure. 

The goal of the dedicated matrix-element-based part is to efficiently generate 
complicated event topologies using native (multi-channel based) phase space 
generation procedures. The strategy is based on the understanding that a case- 
by-case optimisation is in complex cases of phase space topologies preferable to 
an universal algorithm. Given that phase-space is optimised on a case-by-case 
basis, an user-defined pre-selection for the generated regions of the phase-space 
is not implemented. Due to the fact that the 2 — > 4 matrix elements, provided 
by the MADGRAPH/HELAS [6] package, contain full massive treatment of the final 
state particles, there are no explicit divergences present for implemented pro- 
cesses and AcerMC can indeed cover the full (kinematically allowed) phase 
space of the processes at hand. 

The matrix-element-based part uses q:qcd(Q 2 ) and aQED^Q 2 ) couplings and 
mass spectra, as calculated by the supervising event generator, to insure the 
full internal consistency in treatment of the event itself. Optionally, the native 
a QCD(Q 2 ) and olqed{Q 2 ) definitions can also be invoked. 

The generation chain is built from the following steps: 

• The PYTHIA 6 . 2 or HERWIG 6 . 3 interfaces to the library of the structure 
functions PDFLIB 8 . 04 [39] are used to calculate convolution of the partonic 
density. 

• AcerMC modules produce unweighted hard-process events with colour flow 
information and pass them to the supervising generator PYTHIA 6 . 2 or 
HERWIG 6.3 as an external event. 

• The generated events are then further treated within PYTHIA 6 . 2 or HERWIG 
6 . 3 event generators, where the fragmentation and hadronisation proce- 
dures, as well as the initial and final state radiation are added and final 
unweighted events are produced. 

The AcerMC efficiency^ for generating unweighted events, using the imple- 
mentation of the phase-space generation discussed below, is summarised in 
Table 7. A certain (very small) fraction of events is further rejected in the 
showering/fragmentation procedures of the supervising generators. 



6 Note that efficiency is energy scale dependent and phase-space optimisation is 
done individually for each choice. So it might vary for the same process but different 
choices of the energy scale definition. 
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In the following we will briefly describe the key points of the implemented 
AcerMC modules and developed algorithms: matrix element calculations, 
four-fermion phase-space generation, the issue of the s-dependent width for 
resonances, unitary generation of multipheral topology, and finally, the modi- 
fication of the VEGAS algorithm. 
Table 7 

Efficiency for the generation of unweighted events with the default definition of the 
energy scale (see Section 6.6 for details). For generation of the qq,gg — ► Z/7*(— ► 
££)bb andqq,gg — ► Z/j*(^ t€)tt events threshold mm > 60 GeV has been used. For 
generating gg — > [ZjWj^f — >)ttfe& central energy scale was used. The f = e,fj,,T,b. 



Process 


Description 


Internal AcerMC efficiency 


[1] 


gg — > ttbb 


20.2 % 


[2] 


qq — > ttbb 


26.3 % 


[3] 


qq -► 


33.0 % 


[4] 


-► ui)tt 


25.2% 


[5] 


gg^Zh*{^ tt)bb 


33.0 % 


[6] 


qq^ Z/Y(-> U)bb 


29.7 % 


[7] 


gg^Z/j*^ ff,vu)it 


28.2 % 


[8] 


qq^Zh*{^ ff,uu)tt 


34.6 % 


[9] 


gg - (Z/W/7* 


11.2 % 



4-1 The Matrix Element Calculation 

The squared matrix elements of the processes were obtained by using the 
MADGRAPH/HELAS [6] package. They take properly into account the masses and 
helicity contributions of final states particles, incoming quarks are considered 
as massless. The particle masses, charges and coupling values that are passed 
to the code derived with the MADGRAPH package are calculated from functions 
consistent with the ones used in supervising generators (PYTHIA/HERWIG). This 
allows to preserve the internal consistency of the event generation procedure. 
In particular, the (constant) coupling values of a s and «qed were replaced 
with the appropriate running functions that were either taken from the inter- 
faced generators or provided by the AcerMC code according to user settings. 
Slightly modified MADGRAPH/HELAS allowed for obtaining colour flow informa- 
tion of the implemented processes. 

The sets of the MADGRAPH/HELAS coded diagrams, for each of the implemented 
processes, are collected in Appendix A. 
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4-2 The Four Fermion Phase Space Generation 



The four-fermion phase space for the implemented processes was modelled us- 
ing the importance sampling technique based on the procedures implemented 
in the e+e" event generators FERMI SV [35], EXCALIBUR [36] and NEXTCALIBUR 
[37] . For each implemented process a sequence of different kinematic diagrams 
{channels) modelling the expected event topologies was constructed and the 
relative weights between contributions of each sampling channel were sub- 
sequently obtained by using the multi-channel self-optimising approach [38]. 
Eventually, additional smoothing of the phase space was obtained by using a 
modified VEGAS routine to improve the generation efficiency. 

The procedure of multi-channel importance sampling used in the event gen- 
eration can briefly be outlined as follows. An analytically integrable function 
g($), which aims to approximate the peaking behaviour of the differential 
cross-section dependence on various kinematic quantities is introduced into 
the differential cross-section equation as: 

s ($) _ _ _ 

da = s($) rf$ = -W- • g($) d$ = w(<$>) g(<$>) d$, (I) 



where the c?$ denotes the (four-) particle phase space and the s($) summarises 
the matrix element, flux and parton density functions, which all depend on the 
chosen phase space point. The function g(<&) has to be a normalised probability 
density: 

J g($)d$ = 1. (2) 

Since the peaking behaviour of can be very complex due to several pos- 
sible topologies introduced by a large number of contributing Feynman dia- 
grams, the function g($) is composed as a weighted sum of several channels 
<7i($), each adapted to a certain event topology: 

<?($) = (3) 

i 

The values of relative weights ctj are determined from multi-channel self- 
optimisation procedure in order to minimise the variance of the weights 
[38]. The phase space points are than sampled from the function g(<&), first 
by randomly choosing a channel % according to the relative frequencies on and 
then deriving the required four momenta from the chosen gi(<&) using unitaryQ 
algorithms [35]. 

7 Unitary in this context meaning that there is no event rejection in the algorithm. 
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In order to have a closer look at the event generation steps, one first has to 
write down the generic differential cross-section formula: 



|2 



da = ^ f a (x u Q 2 )f b (x 2 , Q 2 )j^^^jdx 1 dx 2 d$, (4) 



where f a ,b( x , Q 2 ) represent the gluon or (anti)quark parton density functions, 
\M.\ 2 the squared matrix element divided by the flux factor 2s and d§ denotes 
the phase space differential. The quantity s = x\ x 2 s is the effective centre-of- 
mass energy, and the sum J2 a ,b runs in case of quark- ant iquark incident partons 
over all possible quark- ant iquark combinations (a, b = u,d,s,c,u,d,s,c). In 
case of gg initial state the sum has only one term with a = b = g. 

Alternatively, in order to use the PDFLIB built-in structure functions xf(x, Q 2 ), 
it is convenient to rewrite the differential cross-section to the form: 

da = J2 x ifa( x i,Q 2 ) x 2 fb (x 2 , Q 2 ) jTrrijl , dy^d$>, (5) 



with the two new variables given by r = x\ ■ x 2 and y = 0.5 log(a;i/x2). The 
matrix element used in the calculation depends explicitly on the four-momenta 
of the incoming and outgoing partons: 

Mab = Mdb(qi,q2,Pl,P2,P3,P4), (6) 



where the Pi represent the four- momenta of the final state particles and qi, q 2 
the four-momenta of the incident partons. Consequently, all the four-momenta 
have to be explicitly generated. 

While generating events with (anti-)quarks in the initial state, an additional 
step is required to pick the flavours of the incoming pair. The selection is again 
done using importance sampling, using as sampling weights: 

A, =x 1 f qi {x 1 ,Q 2 ) -x 2 j- qi [x 2 ,Q 2 ) ■ \Mi\ 2 , (7) 



where the index % runs over all possible quark-antiquark combinations. When 
the matrix element of the hard scattering process is for a given processes 
independent of the incoming quark flavour, it is excluded from the above 
weights for the sake of simplicity (note that the incoming quarks are treated 
as massless so only flavour/(weak-)isospin dependent couplings can introduce 
the flavour dependence of the \Mi\ 2 term). 

Only four flavours (u, d, s, c) of the incoming quarks are considered at the 
moment, the contributions of the incoming b quarks are excluded from cal- 
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culation due to the very high suppression induced by either the structure 
functions and/or CKM matrix suppression. 

In the event generation procedure, a generation channel is thus chosen by 
weighted sampling using a (pre-determined) set of ctj. Next, the values of r 
are sampled from a distribution: 



and y from the distribution 
1 



cosh(z/ y) 



z/~l, (9) 



as also used in PYTHIA standard phase-space algorithm [23]. From the two 
values xi and x 2 , the momenta qi,q2 of incoming particles and the effective 
centre-of-mass energy are derived. In the following step, the four momenta 
of the final state particles Pi are sampled by re-parametrising the general four- 
body phase space: 

d$ = (n <*V(P? - mMpf)) S 4 (qi + g 2 -X»- (10) 
\i=i ) i=i 



in terms of kinematic variables that are expected to exhibit the strongest 
peaking behaviour for a specific channel topology and consequently introduc- 
ing appropriate sampling functions for these variables. 

The modelling of kinematic channels relies heavily on the procedures devel- 
oped in NEXTCALIBUR program [37]; nevertheless, many additions and improve- 
ments were made. Two examples of the extended/added procedures used in 
AcerMC are given below. 

The detailed description of the implementations of four-momenta sampling 
in all existing kinematic channels is omitted for the sake of brevity; the two 
examples below should serve as a representative illustration. 

Subsequently, the four-momenta constructed from the obtained set of kine- 
matic variables are used in the matrix element calculation. Each event is fur- 
ther weighted by the appropriate phase space weight corresponding to the 
importance sampling procedure and calculated using the unitary algorithms; 
for further details on the applied method and unitary algorithms the reader 
is referred to the original papers (e.g. [35,37]). 
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4-2.1 Breit-Wigner Function with s- dependent Width 



In some topologies of the processes involving W or Z° bosons, a bias of 
the matrix element towards large values in the high s^, z region is evident, 
which in turn means that a more accurate description of the tails of s^/z 
distribution is needed. Consequently, the Breit-Wigner sampling function was 
replaced bjf^J 



which is proportional to the (more accurate) Breit-Wigner function with an 
slyy dependent width (W in the above formula denotes either a W or a Z° 
boson). 

In order to implement a unitary algorithm (an algorithm that produces a result 
for every trial, i.e. there is no rejection) of value sampling on the above function 
one first has to calculate the normalisation integral (cumulant) and then its 
inverse function. Introducing a new variable r\ = (sij^ — M^)/(Mw IV) the 
integral of the above function can be expressed as: 

/ BW.m is' w = /{ J^f- . 1 + (12) 



where the upper integral limit is left as a free parameter. The integral thus 
gives a function: 

F(V) = {j^f- ■ atanfa)} + [\ ■ \og( V 2 + 1)}, = F 1 (rj) + F 2 ( V ) (13) 



with F(r] max ) — F(r) m i ri ) defining the normalisation. One of the undesirable 
features is that the function F(r]) does not have a (simple) analytical inverse, 
which is a prerequisite for unitary sampling. Taking a closer look at the two 
above expressions one can quickly spot another undesirable feature, namely 
that the second term in the Equation 12 is an odd function of rj, which after 
the integration gives an even term F2{rj) in rj in Equation 13. In other words 
the second term alone is neither a non-negative function nor does it have 
an unique inverse - one has to deal with a negative probability. A reasonably 
elegant solution to this problem has been developed and implemented here: 

• One samples values of rj by using only the first term of the above expressions 
(the usual Breit-Wigner function). 

8 To our knowledge this implementation is original and done for the first time in 
AcerMC. 
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• One then re-samples the obtained value of rj using the full expression of 
Equation 12: If rj is less than zero the value is mapped to —rj with the 
probability given by Equation 12. 




r| M(W) [ GeV/ c 2 ] 



Fig. 3. Left Comparisons of the two functional terms of Eq. 12 to BW s (t]) given 
by Equation 11. Note that the scaling factor A is chosen in view of making the 
contributions more transparent; it is much too small compared to the real case of 
W ± /Z° bosons. 

Right Comparison of the (normalised) distributions of differential cross-section for 
the process qq — > Wbb (dashed) and sampling functions (solid line) with respect to 
the variables obtained by importance sampling, as described in the text. 

Why this works can quickly be deduced by looking at the Figure 3: At negative 
values of 77 the second term of Equation 12 gives a negative probability in the 
region i] < 0, i.e. using a simple Breit-Wigner (Cauchy) probability function 
too many events are generated in this region. Correspondingly, since the second 
term of Eq. 12 is an odd function, exactly the same fraction (distribution) of 
events is missing in the region rj > 0. By mapping events with rj < over the 
i] = axis one thus solves both problems at the same time. Using the above 
re-sampling procedure the whole approach remains unitary, i.e. no events are 
rejected when there are no limits set on the value of i] or they are symmetric 
|^min| = ^max- In the contrary case, a small fraction of sampling values is 
rejected. 

After some calculation the whole unitary procedure can thus be listed as fol- 
lows: 

• Calculate the kinematic limits 77 min and r) max . 

• Calculate the normalisation factors Ai = Fi(r) max )— iq(?7min), A 2 = F 2 (r] max ) — 
-^K^mm) and A s = A 1 + A 2 ; the term A 2 can actually be negative and thus 
does not represent proper normalisation. 

• Obtain a (pseudo-)random number p\. 

• If pi < A 2 /A s then: 

• Obtain a (pseudo-)random number p 2 ; 
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Construct 77 as: 



X = A 2 - p 2 + F 2 {r] min ), 

which is the inverse of the (normalised) cumulant (F 2 (rj) — F 2 (r] min )) / A 2 . 

■ Note that the condition pi < A 2 /A s can be fulfilled only if A 2 > 0, which 
means that i] max is positive and greater than r] min . 

Conversely, if p\ > A 2 /A s then: 

■ Obtain a (pseudo-)random number p 2 ; 

■ Construct 77 as: 

X = A 1 -p 2 + F 1 (r ]min ), 

,m w t w 

^ = tall ( M 2 - X ) 

1V1 W 

which is the inverse of the (normalised) cumulant {Fi(rf) — i 71 1 (r] min ))/A 1 . 

• If the obtained 77 is less than zero then calculate the normalised probability 
densities: 

P= — f M ™ 1 1 
1 Ax X M W T W l + ri 2 * 

p=— r M ™ 1 1 g 1 

s A s X M W T W 1 + Ty 2 

■ Obtain a (pseudo-)random number p%; 

■ If p 3 > P s /Pi map 77 -> -77. 

• If the new 77 falls outside the kinematic limits [77 min , rj max ] the event is 
rejected. 

• Note also that the last mapping can only occur if the original rj was neg- 
ative, since P s < P± only in the region 77 < 0. 

Calculate the value of s^y using the inverse of 77 definition: 

s* w = (M w T w )- V + M 2 W (14) 
The weight corresponding to the sampled value 77 is exactly: 



(15) 



which is the (normalised) inverse of Equation 11 as requested. 



As it turns out in subsequent generator level studies, this generation procedure 
gives much better agreement with the differential distributions than the usual 
(width independent) Breit-Wigner; an example obtained for the qq — > Wbb 
process is shown in Figure 3. The evident consequence is that the unweighting 
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efficiency is substantially improved due to the reduction of the event weights 
in the high region. 



4-2.2 Unitary Generation of Multipheral Topology 

The procedure implemented in AcerMC for generating events with multi- 
pheral topology is an extension of the three-body phase-space sampling, de- 
veloped in [37] for the phase space generation in NEXTCALIBUR program (kine- 
matic channel MULTI1, Appendix A). In order to describe the features imple- 
mented in the extension of that procedure, let us recapitulate first its basic 
principles: The procedure used in MULTI1 starts with splitting the four-body 
phase space into a three-body times a two-body decay: 

J rf$ 4 = J ds u h(pi,P2,P34) h(p3,P4), (16) 



where the integral 7 3 represents the three-body and I2 the two-body phase 
space, the labels Pi, % — 1,4 represent the four-momenta of the produced par- 
ticles with p 34 = P3+P4 and s 34 = pl 4 . The key issue thus becomes parametri- 
sation of the (more difficult) three body phase-space integral I 3 (Pi, P 2 , P3): 

h(Pi, P2, Pb) = I f[ d * P ^P- - Mf)e(E i ). (17) 
J i=i 



It should be stressed that the masses Mj in the above on-shell constraint do 
not necessarily correspond to physical masses of the elementary particles, can 
as well be values derived by a pre-sampling procedure (e.g. Mf = s 34 in the 
above example, M 3 representing the off-shell mass of a gluon decaying into p 3 
and p^. Defining the reduced energies x and y as: 

Ei = ^x, E 2 = ^y, (18) 



where the s is the square of centre-of-mass energy of the system, and setting 
Hi = Mf /s one can re-parametrise the integral h{P\, P2, P3) as: 

X- 

/ 3 (Pi,P2,P 3 ) = ^ / dx j 'dy<m l( m25[F(x : y)}. (19) 



The integral is now parametrised with the energies and solid angles of two 
particles with a constraint in form of a S[F(x,y)] due to the presence of the 
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third particle. The integration limits on x can be easily derived from the Dalitz 
constraints, with x_ given by: 



X- = 



(1 - (v^+ y/JI^) 2 . 



(20) 



The delta-function 5[F(x, y)\ is actually a constraint on the value of the cosine 
of the angle between the particles 1 and 2: 



the full expression for C(x,y) is given in [37]. A short calculation shows that, 
given values x G [2^///!, x_] and C12 € [— 1, 1], the constraint of Equation 21 
gives either two, one or no solutions for y, i.e. the procedure of event generation 
can involve some event rejection and is thus not necessarily unitary. 



Fig. 4. A diagrammatic representation of a three-body decay as modelled by the 
procedure described in the text. 

In the original MULTI1 channel implementation, the event generation uses the 
sampling sequence: 

• Sample the value of x. 

• Calculate the value of y assuming cyi = —1 (which is favoured by multi- 
pheral singularity). 

• Use the (x,y) to sample the four values cos^i^), 0i,2- 

• From the angles calculate c±2 = cos(#i2) and hence determine the possible 
values of y using Equation 21. If there is no solution for y, the event is 



As it turns out, the rejection rate, very low in cases of processed where it 
was applied originally [37], increases dramatically for the QCD process where 
heavy quarks are involved. Therefore, the approach was extended to an unitary 
one (albeit at the cost of computation time). 



F(x,y) = C(x,y) - c 12 , 



(21) 




rejected. 
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The developed extension is based on a short calculation, which shows the 
existence of a turning point (maximum) at i/t = Vt{x) m the Ci2(y) = C(x,y) 
function, i.e. there is a value of x where Cuiur) reaches its maximal value. At 
given values of x and Ci 2 , the number of solutions for y thus depends on the 
location and value of the extremum ct = cuiyr)- Basing on these deductions, 
a new approach was developed: 

• Sample the value of x from a distribution: 

' (22) 



(1 -xf 

• Sample the angles <f>i G [0, 27r] from a flat distribution and 6 X from: 

1 „ 2E l E 9l - Si- S gi + S 3 

77TT a l ~~ 1 - 1 1 ~, > 

ai-cos(0i) 2|Pi||P Pl | 

where the index g\ represents the incoming particle that splits into the t- 
channel exchange and the Sj terms correspond to Sj = M 4 2 from Equation 
17 (remember that E\ = \fsj1 ■ x). 

• Calculate the values of kinematic limits on [y m m{x) , y ma ^(x)] and the values 
yrix). If the yr is within kinematic limits [y m i n (x) , y mSuX (x)], calculate the 
maximum limit Ct and hence the limits on 2 and 2 . 

• Calculate the value of y assuming c i2 — — 1 (which is always allowed). 

• Sample the value of 9 2 within the calculated limits from a distribution anal- 
ogous to Equation 23 using y. 

• Get the value of 02 from a flat distribution within the allowed limits. 

• From the generated angles calculate c 12 — cos(#i 2 ) and hence determine 
the possible values of y using Equation 21. If the calculations were done 
correctly, there should always be at least one solution. 

The modified procedure is thus an unitary one, giving the four-vectors Pi, P 2 , P 3 . 



Going back to the initial issue of the four body phase space (Eq. 16), the whole 
sampling procedure for such kinematic channels is as follows: 



The squared invariant mass S34 is sampled from: 
1 



(24) 



The value of \x depends on the channel implementation (virtual gluon, a 

t-channel object, masses of the particles 3 and 4). 

The described above three-body decay procedure is implemented. 

The virtual P 34 object is decayed according to the actual channel (a t- 

channel angular dependence, isotropic decay of a virtual gluon). 
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4-2.3 Modified VEGAS Algorithm 



Using the described multi-channel approach, the total generation (unweight- 
ing) efficiency amounts to about 3 — 10% depending on the complexity of the 
chosen process. In order to further improve the efficiency, a set of modified 
VEGAS [5] routines was used as a (pseudo-) random number generator for sam- 
pling the peaking quantities in each kinematic channel. The conversion into a 
(pseudo-)random number generator consisted of re-writing the calling routines 
so that instead of passing the analysed function to VEGAS for sampling and in- 
tegration, VEGAS calls produce only (weighted) random numbers in the region 
[0, 1] and the corresponding sampling weight, while the VEGAS grid training is 
done using a separate set of calls. 

After training all the sampling grids (of dimensions 4-7, depending on the 
kinematic channel), the generation efficiency increased to the order of 6 — 14%. 
The motivation for this approach was that in unitary algorithms only a very 
finite set of simple sampling functions is available, since the functions have 
to have simple analytic integrals for which an inverse function also exists. 
Consequently, the non-trivial kinematic distributions can not be adequately 
described by simple functions at hand in the whole sampling domain (e.g. the r 
distribution, c.f. Figure 5) and some additional smoothing might be welcome. 
In addition, the random numbers distributions should, due to the applied 
importance sampling, have a reasonably flat behaviour to be approached by 
an adaptive algorithm such as VEGAS0. 

The further modification of VEGAS, beside adapting it to function as a (pseudo- 
random number generator instead of the usual integrator, was based on the 
discussions [41,42] that in case of event generation, i.e. unweighting of events 
to the weight one, reducing the maximal value of event weights is in principle 
of higher importance than achieving the minimal weight variance. Since the 
VEGAS algorithm was developed with the latter scope, some modification of the 
algorithm was necessary. As it turned out, the modification was fairly easy to 
implement: Instead of the usual cumulants: 

</>cell =£wti, (25) 
cell 



according to the size of which VEGAS decides to split its cells, the values: 

< F > ce ll = Accu ■ wt™ x - J2 wti, (26) 

cell 

9 At this point also a disadvantage of using the adaptive algorithms of the VEGAS 
type should be stressed, namely that these are burdened with the need of training 
them on usually very large samples of events before committing them to event 
generation. 
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were collected and used as the splitting criterion. The above value (called loss 
integral in [42]) is basically a measure of the deviation between the maximal 
weight sampled in the given cell wt™ff and the average weight in the cell 
< wt ce u >= (Xccii wti) / A ce ii (the quantity A ce n denoting the cell width, i.e. 
the integration range). Re- writing the above expression as: 

<F> cell =(A ceU .wCT-{l-^^} (27) 



clearly indicates that the value < F > ce u is actually a measure of the genera- 
tion inefficiency in the cell, since the term in the curly brackets is equivalent to 
one minus the generation efficiency < wt ce n > /wt™ff. In addition, the ineffi- 
ciency is weighted with the crude/ maximal estimation of the function integral 
over the cell A cell • wt™Jf and cells with the highest < F > ccll are split. 

This method is of relevance because the VEGAS cells are actually projections 
of the whole phase space on the (chosen) side axes, i.e. VEGAS cannot isolate 
a maximal weight in a certain point in phase-space and build a cell around it, 
which in principle would be an ideal solution. An implementation with this 
scope in view has been made in FOAM [42], nevertheless we have not found it 
competitive with respect to the modified VEGAS for the given application. 

The thus modified ac-VEGAS algorithm further increased the unweighting ef- 
ficiency for almost a factor of two. 



-2 -2 
xlO xlO 




Fig. 5. Comparison between the sampling distribution for the r = s/s £ [r m i n , 1] 
variable in gg— > ttbb process before and after the application of modified ac-VEGAS 
[5] smoothing procedure (light gray histogram). The generated (normalised) differ- 
ential cross-section is also drawn (black histogram, labelled Generated ). 



One of the sampling distributions is shown in Fig. 5 as a gray histogram 
(marked channel) and the actual {generated) differential cross-section depen- 
dence is drawn in black. In the first figure, the random variable used for 
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sampling values from l/r M distribution was drawn from a flat probability in 
the interval [0, 1]; in the second plot the ac-VEGAS algorithm was used to give 
an optimal grid for sampling the random variables needed for parameter gen- 
eration (the grid is trained for each kinematic channel separately, the sum of 
all channels is shown in the plot). The improvement is evident; one has to 
stress that the use of ac-VEGAS algorithm to generate the values of r directly 
would be much less efficient since VEGAS gives a grid of 50 bins/dimension, 
which would give a very crude description of the r distribution compared to 
the one at hand. 




Multi-Chan only 
Multi-ChaneVEGAS 
Multi-Chan0ac-VEGAS 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Wtx5 • 10 



Fig. 6. The distribution of event weights using only the Multi-Channel approach 
(dotted histogram) and after application of VEGAS (dashed histogram) and ac-VEGAS 
(full histogram) algorithms in the gg — ► (Z° —*)llbb process. 



Observing the distributions of the event weights before and after the inclusion 
of the modified ac-VEGAS algorithm (Fig. 6) it is evident that ac-VEGAS quite 
efficiently clusters the weights at lower values. Note that the principal effect 
of original VEGAS is indeed to cluster event weights in a narrow region, never- 
theless a tail towards the high-weight region remains. On the other hand, the 
ac-VEGAS efficiently reduces the tail in the high weight region; only a few of 
the event weights still retain their large values, thus reducing the generation 
efficiency. Given the difference in distributions, the observed increase of the 
generation efficiency seems relatively modest. To better understand this result 
one should consider that the formula for the MC generation efficiency is given 
by: 



< wt > 

Wtmav 



(28) 



where < wt > is the average weight of the sample and equals the total event 
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cross-section, while wt max represents the maximum event weight in the applied 
generation procedure and is determined through a pre-sampling run with a 
high statistic. Since the average weight < wt > equals the total cross-section 
of the process, it remains (necessarily) unchanged after the application of the 
VEGAS refining; consequently the change of efficiency results in the reduction 
of the maximum weight wt max by approximately a factor two, which is from 
technical point of view quite an achievement. 

A further step to profit from the clustering of weights induced by ac-VEGAS is 
to adopt a re-definition of the MC generation efficiency as proposed by [40,42]. 
In this approach, the alternative definition of wt max is: For a given precision 
level a « 1, the wt max is determined from the total weight distribution in 
such a way that the contribution of the events exceeding this value to the total 
weight sum (i.e. cross-section integral) equals a. Such a quantity is referred 
to as wt max and the efficiency expression becomes: 

* = (29) 
wt a 

max 



The argument presented in [40,42] seems to be quite reasonable since the true 
event weight is in any case only estimated from a finite sample of events and 
the new definition simply takes into account a certain level of accuracy in the 
maximum weight determination. In addition, certain very weak singularities 
that might exist in the simulated process and might occasionally result in a 
very high event weight are automatically taken into account. The use of new 
wt max consequently results in a generation efficiency of about e > 20% for 
all the implemented processes, which is a significant improvement in terms of 
time needed for MC generation. 



4-2.4 Colour Flow Information 

Before the generated events are passed to PYTHIA/HERWIG to complete the 
event generation, additional information on the colour flow/connection of the 
event has to be defined. Below we discuss the implemented method of the 
colour flow determination on the example of two processes, gg — > tibb and 
qq — > tibb. 

For the process gg — > tibb six colour flow configurations are possible, as shown 
in Figure 7. With 36 Feynman diagrams contributing to the process and at 
least half of them participating in two or more colour flow configurations, cal- 
culations by hand would prove to be very tedious. Consequently, a slightly 
modified colour matrix summation procedure from MADGRAPH [6] was used to 
determine the colour flow combinations of the diagrams and the corresponding 
colour factors. The thus derived squared matrix elements for separate colour 
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flow combinations |-Mfl ow | 2 were used as sampling weights on an event-by- 
event basis to decide on a colour flow configuration of the event before passing 
it on to PYTHIA/HERWIG for showering and fragmentation. The procedure was 
verified to give identical results regarding the colour flow combinations and 
corresponding colour factors when applied to the processes published in [43]. 
As one can see this approach neglects the interference terms between the dis- 
tinct colour-ordered amplitudes and is indeed exact only in the Nc — > oo 
limitf 1 '] [44,45]. Since there is no a priori rule of how to split the interfer- 
ence terms between the colour ordered amplitudes this approach is generally 
deemed to be the best one can do; recent developments in this field [46] how- 
ever suggest additional improvements to the method that indeed might be 
incorporated into later versions of AcerMC. 




Fig. 7. A diagrammatic representation of the six colour flow configurations in the 
process gg — > tibb. Certain colour combinations, leading for example to colourless 
(intermediate) gluons, are not allowed. 

The colour flow configuration in the qq — > tibb channel is much simpler since 
only two colour flow topologies exist (Fig. 8); the choice between the two has 
been solved in a manner identical to the one for the gg — > tibb process, as 
described above. 




Fig. 8. A diagrammatic representation of the two colour flow configurations in the 
process qq — > tibb. 

10 The matrix elements used in the cross-section calculation and event generation 
are of course complete and do not employ any approximation. 
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Some specifications of the implemented matrix-element-based processes: num- 
ber of Feynman diagrams, channels used in the phase-space generation and 
colour flow configurations are collected in Table 8. 

Table 8 

Some details on matrix-element-based process implementation in AcerMC library. 
In case of qq initial state the number of Feynman diagrams corresponds to one 
flavour combination. The f = e,/i,T,b. 



Process id 


Process specification 


Feyn. diagrams 


Channels 


Colour flows 


1 


gg — > ttbb 


36 


11 


6 


2 


qq — > ttbb 


7 


4 


2 


3 


qq^W{^ lv)bl 


2 


2 


1 


4 


qq -» W(-> tv)tt 


2 


2 


1 


5 


gg ^Zh*(^U)bb 


16 


6 


2 


6 


qq^ Z/7*(-» ££)&& 


8 


5 


1 


7 


gg^Z/>f{-+ ff,uu)ti 


16 


6 


2 


8 


qq^Z/ 7 *(^ ff,vu)tt 


8 


5 


1 


9 


gg -> (Z/VT-/7* ^)tt"66 


72 


21 


12 



^.5 The ckqed and ot s calculations 

Native functions of running ojqed(Q 2 ) and a s (Q 2 ) have been implemented 
inside AcerMC with the main objective of providing a means to keep the 
(total) cross-sections of the processes unchanged when interfacing with the 
two supervising generators, since the implementations of the two functions in 
PYTHIA and HERWIG differ to some extent. Especially the a s (Q 2 ) is subject 
to experimental and theoretical uncertainties, however obtaining a different 
cross-sections for the same AcerMC process due to different interface, could 
be regarded (at least to some extent) as an inconsistency^]. 

• «qed is implemented in AcerMC using the formulae given in [47] and is in 
complete accordance with the implementations in PYTHIA and HERWIG apart 
from the updated hadronic component published recently by Burkhardt et. 
al. [48]. As one can see in Figure 9, the latter minimally lowers the «qed 
values. 

11 The values will still differ by a small amount in processes containing W bosons 
(processes 3,4) due to different values of the CKM matrix in the two supervising 
generators. 
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Fig. 9. Comparison between the chqedIQ 2 ) (Left) and a s (Q 2 ) (TiightJ implementa- 
tions in AcerMC, PYTHIA and HERWIG. For a s (Q 2 ) calculations with different loop 
orders (L) are given where applicable. 

• a s has one and three loop implementations in AcerMC following the cal- 
culations of W. J. Marciano [49] and using A^{? transformations for flavour 
threshold matches. The three loop version gives good agreement with the 
HERWIG implementation (both functions have been set to the same A^ -5 ^ 
value) as one can see in Figure 9. The PYTHIA two loop implementation 
deviates somewhat from the latter two; the kinks observed in the plot are 
due to approximate A^{? transformations at flavour thresholds, which are 
exact to one loop only. 

Although the AcerMC and PYTHIA one loop implementations are identical in 
form the resulting values differ by a small amount because the default PYTHIA 
implementation reads the A^ -4 "* value from PDFLIB instead of the A^ _5 ' ) 

one used by AcerMC and HERWIG; the difference thus occurs due to A^{? 
propagation at flavour thresholds. 
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5 Structure of the package 



The AcerMC package consist of a library of the matrix-element-based genera- 
tors for selected processes, interfaces to the PYTHIA 6 . 2 and HERWIG 6 . 3 gen- 
erators, sets of data files and two main programs: demojiw.f and demo_py.f. 
Provided makefiles allow to build the executables with either of these genera- 
tors as the supervising generator, demo_hw.exe and demo_py.exe. 

5.1 Main event loop and interface to PYTHIA/HERWIG 

The main event loop is coded in the demojiw.f or demo_py.f files, where 
the opening/closing of the input/output files, reading of the data-cards and 
event-loop execution is performed. Main event loop consists only of calls to 
the acermc_py or acermc_hw subroutines, with parameter MODE = -1, 0, 1 
respectively set for initialisation, generation and finalisation of the event loop. 
The call to acermc_xx activates respective procedures of the supervising gen- 
erator, which in turn activates the acevtgen procedure steering the native 
AcerMC generation of the matrix element event. Fig. 10 illustrates this call- 
ing sequence in some details. 

As one can deduce from the diagram in Fig. 10, certain functions called by 
AcerMC, as e.g. pseudo-random number generator acr are re-routed through 
the interfaces to the linked supervising generator, depending on the choice 
at compilation time (e.g. acr function giving (pseudo-)random numbers is 
linked to either pyr or hwrgen as shown in the plot), providing the internal 
consistency of the package. The generated event is rewritten to the format 
required by the supervising generator by means of the acdump_xx routines. 

The pythia_ac.f and herwig_ac.f files contain sets of re-routing/interface 
functions, specialised for the respective supervising generator. The main li- 
brary of AcerMC is well screened from dependencies on the supervising gen- 
erator, all dependencies are hidden in herwig_ac.f and pythia_ac.f respec- 
tively. 

The PYTHIA and HERWIG libraries remain essentially untouched^ 2 ], without 
introducing any dependencies on the AcerMC code. The input cards are 
common for both interfaces. 



38 

For specification of exceptions see Section 6.4 and 6.5. 



demo_py . exe 



acermc _py(-lh 
acermc _py(0)- t 
acermc _py( 1) -] 



Input cards 
run. card 
acermc.card 



demo_hw . exe 



acermc_hw( 
acermc_hw(0) r 
acermc _h\v( 1 ) p 



Pythia 6.2 



hwepro i(— 



hwefin 



pyr 
pypdfu 
pyalem 
pyalps 



pythia_ac 



Maxima^ 
weights 



ac-VEGAS 
grids 



5 



acr 
acpdf 
ac_alpem 
ac_alps 



acdump_xx 



Chanel 
weights 
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acevtgen 



herwig_ac 
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▼ i 



acevent_xx 



RcerMC 1 . 




HERWIG 6.3 



hwigin 
hwuinc 
hweini 



hwrgen 
hwsfun 
hwuaem 
hwualf 




Fig. 10. The calling sequence of the main event generation routine acevent_xx. 
The routine is called either through demo_py — > acermc_py sequence when inter- 
facing the PYTHIA 6.2 generator or demoJiw — > acermcJiw sequence when the 
HERWIG 6 . 3 is linked. The structure of the interface subroutines and relations with 
the corresponding ones from supervising generators and/or external libraries is also 
evident. 



5.2 Structure of the AcerMC matrix- element and phase-space code 



The AcerMC core code performs the generation of a matrix-element-based 
event. Fig. 11 illustrates an example of the calling sequence for generating 
gg — > tibb event. The steering subroutine is called acevent_xx, where xx 
denotes an unique label corresponding to the process at handf 13 ] (in case of 
gg — > tibb we have xx=tt). To stress again, this subroutine calls only a se- 
quence of the native AcerMC subroutines, any call to the supervising gen- 
erator goes via the respective interface function/subroutine. A more detailed 
representation of calling sequence is shown in the Figure 11. 



3 The label does actually not have any relation with the process at hand, it is just 
an unique two character choice. 
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Fig. 11. The event generation sequence controlled by acevent_tt subroutine. Phase 
space generation is sequenced by calling first the acshat_tt routine to obtain the 
incoming gluon momenta and next the phase space generation via the ackin_sel_tt 
routine which re-routes the generation to a certain kinematic channel ackin_ch. 
The latter routine handles the possible momenta permutations and calls the explicit 
four-momenta generation (and PS weight calculation) in the routines acwtgen_tt. 
These (channel-specific) routines are constructed from common building blocks listed 
in the next two columns. The acevent_tt routine also initialises MADGRAPH/HELAS 
package and retrieves the matrix element values. All the generated four-momenta, 
as well as the event weight are finally passed back to the supervising generator via 
the acdump_tt call. 

Code for the phase space generation is grouped together in the subdirectories, 
one per subprocess, e. g. code for generating gg — > ttbb event is in subdirectory 
01_gg_ttbb. Code for matrix element calculations is grouped together for all 
processes in subdirectory matel. Code with different utility subroutines, e.g. 
kinematic transformations used by all subprocesses, is in the subdirectory 
common. Subdirectory interface contains code with interfaces to supervising 
generators, finally subdirectory include contains all include files. The overall 
view on the structure of the AcerMC directories is shown in Fig. 12. 



The core code builds one library libacermc.a. 



5. 3 Data files for the phase-space optimisation 



The AcerMC matrix-element-based generators are very highly optimised, us- 
ing multi-channel optimisation and additional improvement with the ac-VEGAS 
grid. The generation modules require three kinds of the input data to perform 
the generation of unweighted events: 
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utils 



Fig. 12. The structure of the AcerMC directories. 

• A file containing the list of the values of relative channel weights obtained by 
the multi-channel optimisation, defaults being stored in acermc_src/ include. 

• A file containing the pre-trained ac-VEGAS grid, the pre-trained (default) 
ones located in acermc_dat/grids. 

• A file containing the maximum weight wt max , a-cutoff maximum weight 
wt^, ax and the 100 events with the highest weights, the default ones being 
provided in acermc_dat/maxims. 

In case of changing the default running conditions, like parton density func- 
tions or centre-of-mass energy, the user should repeat the process of prepara- 
tion of the listed data files containing the inputs for the phase-space generator 
modules in order to preserve the initial event generation efficiency. 

The reading sequence of data files inside AcerMC is shown in Figure 13. 
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Fig. 13. The reading sequence of the input files and the performing subroutines in 
the AcerMC code. 



Pre-trained data sets are obtained using ^Js = 14 TeV, PYTHIA default a s (Q 2 ) 
and aQED(Q 2 ) and CTEQ5L (parametrised) parton density function set and 
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are provided for each implemented process^]. For these, the relative channel 
weights are stored in the INCLUDE files in acermc_src/ include/ chanwt_xx . inc 
where xx denotes the process id (c.f. Table 7); the default /pre-trained ac-VEGAS 
grids are listed in the directory acermc_dat/grids/vscalA_xxYYY. veg, where 
A denotes the scale choice of the process xx and YYY denotes the cutoff value 
of the m^o/ 7 * for the AcerMC processes xx = 05 — > 08. The files containing 
the maximal weights wt max and wt^, ax as well as the 100 events with the high- 
est weights are stored in the directory acermc_dat/grids/vtmaxA_xxYYY.dat, 
following the same labelling convention. Both the trained ac-VEGAS grids and 
the weight files were obtained from test runs with at least 2 ■ 10 6 weighted 
events being generated. 

The number of required input files might at first look seem large, considering 
that many event generators do not require any input files for operation; the 
difference is not in so much in the complexity of the phase space generation 
as in the fact that many event generators require a warming run instead, i.e. 
before the generation of unweighted events is performed a certain number 
of weighted events (typically of the order of 10 4 ) is generated in order to 
obtain the relative multi-channel weights (in case multi-channel phase space 
generation is used) and/or the optimised VEGAS grid and/or an estimate of 
the maximal weight. Such an approach can have an advantage when event 
generation is very fast and the phase space regions with the highest weights 
are well known (as done for the 2^2 processes in PYTHIA); on the other 
hand, when the phase space topology of the process is more complex and the 
event generation is comparatively slow, generating a relatively small number 
of e.g. 10 4 weighted events every time a generator is started can become CPU 
wasteful and/or inaccurate in terms of maximum weight estimation. 

Reasonably accurate estimation of the latter is namely crucial for correct event 
unweighting; event generators using warming-up method for maximal weight 
search often find still higher weights during the production run and reset the 
maximal weight accordingly. In this case however, statistically correct ap- 
proach would be to reject all events generated beforehand and start the event 
generation anew, which is almost never implemented due to the CPU consump- 
tion and the possibility of hitting a weak singularity (the same argument leads 
to the definition of the wt^ ax , c.f. Section 4.2.3). With a small pre-sampled set 
the generator can however badly under-estimate the maximum weight and a 
large number of events can be accepted with a too-high probability. The only 
hope of obtaining correct results is in such cases that the weight plateau will 
be hit sufficiently early in the event generation process. Consequently, such 
approach can be very dangerous when generating small numbers of event sp 5 ]. 



These can also be used for a series of other settings, see Section 6.3 for details 
15 Small being a somewhat relative quantifier, since the size of an representative 
sample should depend on the phase space dimension, i.e. the number of particles 
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In contrast to the warming-up approach the AcerMC we decided that using 
separate training runs with large numbers of weighted events to obtain the 
optimised grids and maximum weight estimates are preferable, in case user 
wants to produce data sets for non-defaults setting, this can easily be done by- 
configuring the switches in the acermc.card (see Section 6.2). 



in the final state; with e.g. 4 particles in the final state, 10 5 events can still be 
considered a relatively small statistics. 
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6 How to use the package 

There are two steering input files: run. card and acermc.card which share 
a common format for both executables. The run. card (see App. B.l) pro- 
vides switches for modifying: generated process, number of events, parton 
density functions, predefined option for hadronisation/fragmentation in the 
supervising generator, random number, etc.. The acermc.card (see App. B.2) 
provides switches for modifying more specialised settings for the AcerMC li- 
brary itself. Once the user decides on a setup for the generated process, only 
run. card is very likely to be modified for the job submission. Both input 
files are read by AcerMC executables through the CERNLIB FFREAD routines, 
some commands given in the input files (e.g. LIST entry, see Appendix B) are 
internal FFREAD commands which should be disregarded by the user. 

The same executables can also be used for running standard PYTHIA 6 . 2 and 
HERWIG 6.3 processes. The example how to require such process is provided 
as well, in demojiw.f and demo_py.f respectively. If the user requires that 
the AcerMC library is not used, the ttH production will be generated with 
demo_py.exe and HERWIG 6.3 implementation of the Zbb production will be 
generated with demo_hw . exe. In this case only the run . card file will be read, so 
in case the user requires different processes and/or settings of the supervising 
generators the user has to implement her/his steering there or create another 
xxx.card file, together with the corresponding code added to the demo_xx.f. 



6. 1 Steering switches of the overall run 

The overall run is controlled by the switches read from the run . card file (see 
also App. B.l). Some of these general switches are also passed to the AcerMC 
library. 

• CMS : Sets the centre-of-mass energy in GeV. 

• ACER : Specifies if the internal AcerMC process will be used 
ACER=0 - use process from PYTHIA/HERWIG 

ACER=1 - use internal AcerMC process 

• PROCESS : Sets process id 

• HAD : Sets predefined option for QCD ISR/FSR and hadronisation 
HAD=0 - only hard process 

HAD=1 - only ISR (works for PYTHIA interface only) 
HAD=2 - only ISR and FSR 
HAD=3 - full treatment 
HAD=4 - only FSR 

HAD=5 - only FSR and hadronisation 
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• PDFLIB804 NGROUP 

Sets the value of the PDFLIB804 NGROUP parton density function choice. 

• PDFLIB804 NSET 

Sets the value of the PDFLIB804 NSET parton density function choice. 

• RSEED : Choose the random seed for (pseudo-)random generator initialisa- 
tion 

• NEVENT : Required number of generated events 
6.2 Steering switches of the AcerMC processes 

The AcerMC processes are controlled by values set in a simple arrays speci- 
fied in acermc_src/include/AcerMC . inc: 

C CROSS-TALK PARAMETERS 
DOUBLE PRECISION ACSET 
INTEGER IACPROC 

C0MM0N/ACPAR1/ACSET (200) , IACPROC (200) 

C PARTICLE PROPERTIES 

DOUBLE PRECISION ACCHG , ACMAS , ACCKM 

C0MM0N/ACPAR2/ACCHG(50,4) ,ACMAS(50,4) ,ACCKM(4,4) 

C ROUTINE I/O 

INTEGER LACSTD,LACIO 

C0MM0N/ACPAR3/LACSTD , LACIO 

The IACPROC array activates the process IPR0C=PR0CESS (read from run . card 
file) by setting IACPROC (IPR0C) = 1. 

The list of currently implemented processes in AcerMC can be found in Table 
7. When running in the generation mode with ACER=0 full list of processes 
implemented in either PYTHIA 6 . 2 or HERWIG 6 . 3 can be activated, however 
the mechanism for passing information about process id to either of these 
generators has to be coded by user individually in demo_xx.f. 

The main control switches reside in the array ACSET. The COMMON block ACPAR2 
contains the particle charges, masses and decay widths as well as the CKM 
matrix using the PYTHIA convention. The values are filled by the interface 
routines to be equal to the PYTHIA/HERWIG internal values in order to pre- 
serve consistency within the generation stream. In case the user wants to 
change some of the particle properties this should be done through the native 
PYTHIA/HERWIG switches; AcerMC will copy them and use the new values. 

The COMMON block ACPAR3 contains the two logical I/O unit numbers used 
by AcerMC. The LACSTD value determines the output unit of the AcerMC 
messages and the LACIO unit is used for reading/writing the AcerMC data 
files. 
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The main control switches which reside in the array ACSET (see also App. B.2): 

• ACSET (1) : Sets the centre-of-mass energy in GeV. 

• ACSET (2) : Scale of the hard process 

Choose the Q 2 scale for the active AcerMC process. The implemented 
values differ for various processes, the currently implemented settings are 
specified in Section 6.6. 

• ACSET (3) : Fermion code 

The flavour of the final state fermions produced in W ± , Z°/Y — > // decays 

of AcerMC processes 3^8. The PYTHIA/PDG naming convention is used: 

ACSET(3)=11 - W -> ev e ; Z°/Y -> e+e" 

ACSET(3)=12 - Z°/Y -> v e v e) 

ACSET (3) =13 - W -> fiv^ Z°/Y -> 

ACSET (3) =14 - Z°/Y -> 

ACSET (3) =15 - W -> tv t ; Z°/Y -> t+t~ 

ACSET (3) =16 - Z°/7* -*_i/ T z/ T 

ACSET (3) =5 - Z°/7* -> ^ 

At present the ACSET (3) =5 and ACSET (3) =12, 14, 16 are implemented only 
for processes 7^8. 

• ACSET (4) : Z°/~f* propagator 

Use full Z°/7* propagator instead of the pure Z° propagator in matrix ele- 
ment calculation for the AcerMC processes 5^8. The switch is provided 
since in some of the analyses the 7* contribution is of relevance in the se- 
lected mass windows; for the analyses selecting the mass window around 
the Z° peak this contribution can safely be neglected. 
ACSET (4) =0 - only Z° propagator. 
ACSET (4) =1 - full Z°/Y propagator. 

• ACSET (5) : m^o/ 7 * mass cut 

Cutoff value on the invariant mass m z o/-y* in GeV when ACSET (4) =1. Note 
that the provided data files exist only for values of ACSET (5) =2,5,10,15,30,60, 
120,270,300 and 500 GeV which should satisfy most user requirements for 
the analyses foreseen at LHC. In different value is set the user has 

also to provide the user data files for the run. 

• ACSET (6) : Sets the value of the PDFLIB804 NGR0UP parton density function 
choice. 

• ACSET (7) : Sets the value of the PDFLIB804 NSET parton density function 
choice. 

• ACSET (8): The implementation of a s (Q 2 ) 

Selects the implementation of a s (Q 2 ) to be used in the matrix element 
calculation: 

ACSET (8) =0 - Use the a s (Q 2 ) as provided by the supervising generator/ 
ACSET (8) =1 - Use the a s (Q 2 ) (one loop) provided by the AcerMC; this 
option gives a s (Q 2 ) values equal to the default PYTHIA implementation. 
ACSET(8)=2 - Use the a s (Q 2 ) (three loop) provided by the AcerMC. 
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ACSETO): A ( ™/ s ' value 

Sets the A^T 5 ' ) value to be used in the a s (Q 2 ) calculations in case the Ac- 
erMC native implementation (ACSET(8)=l) is used. 

ACSET(9)=-1 - The A^ =5) value is taken from the PDFLIB804 for the se- 
lected parton density function set. 
ACSET(9)>0 - The provided value is taken. 
ACSET(IO): The implementation of a QED (Q 2 ) 

Selects the implementation of «qed(Q 2 ) to be used in the matrix element 
calculation: 

ACSET(10)=0 - Use the «qed(Q 2 ) as provided by the supervising generator. 
ACSET(10)=1 - Use the a QED (Q 2 ) implemented in the AcerMC. 
ACSET(ll): «qed(0) value 

Specifies the value of q;qed(0) for AcerMC «qed(Q 2 ) calculation. 
ACSET(11)=-1 - The «qed(0) value is set to a QED (0) = 0.0072993. 
ACSET(11)>0 - The provided value is taken. 
ACSET(12): Decay mode of the produced ti pair 

Sets the decay mode of the W boson pair from the ti final state in the 
AcerMC processes 1,2,4,7,8 and 9. For ACSET(12)>0 the combinatoric 
value of the a x BR is recalculated and printed in the output. This switch 
was implemented since the supervising generators (PYTHIA/HERWIG) do not 
allow for forcing specific decays of the top quark pairs generated by external 
processes. This switch imposes a modification of the decay tables of the 
supervising generators on an event by event basis. 

ACSET(12)=0 - both W bosons decay according to PYTHIA/HERWIG switches. 

ACSET(12)=1 - W 1 -> ev e and W 2 -> qq. 

ACSET(12)=2 - Wx -> pv^ and W 2 -»• qq. 

ACSET(12)=3 - Wi -> rv T and W 2 -> qq. 

ACSET(12)=4 - Wi -> ev^iWp and W 2 -> qq. 

ACSET(12)=5 - Wi -> ev^iWp and W 2 -> qq. 

The setting ACSET(12)=5 works for PR0CESS=4 only and implies leptonic 
decay of the VU-boson with the same charge as the one of the primary W 
boson produced in the hard process. Folowing configurations are posible: 

qq -> W + ti -> (W + ^)L + v L (Wf b (W 2 ->)q'q' b, 

or: 

qq -> W'tt -> (W~ ->)L-D L {W? ^)q'q b (W 2 ->)Z - Pj b, 

where L* 1 is the lepton from the primary W decay (controlled by ACSET(3) 
switch) and is either an e ± or ^ as for ACSET(12)=4. Since the charge of 
the semi-leptonic decaying W is correlated with the charge of the primary 
W boson, the a x BR is consequently a factor two smaller than the one for 
ACSET(12)=4. 
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• ACSET(50): AcerMC training mode 

The switch controls the mode in which AcerMC is run: 
ACSET(50)=0 - production run, generate unweighted events. 
ACSET(50)=1 - perform multi-channel optimisation and output the user file 
with channel weights. 

ACSET(50)=2 - perform ac-VEGAS grid training and output the user file with 
trained ac-VEGAS grid. 

ACSET(50)=3 - perform ac-VEGAS grid training as in ACSET(50)=2 but do 
this by updating a provided grid. 

• ACSET(51): Required number of generated events NEVENT 

In case the switch ACSET(50) is set to the non-zero value (i.e. in one of 
the training modes) the ACSET(51) entry is used and defines the number of 
(weighted) events that will be generated; this information is necessary for 
the learning algorithms to decide on steps in the learning sequence. 

• ACSET(52): User data files 

Use the data files provided by user: 

ACSET(52)=0 - no, use native (default) AcerMC data files. 

ACSET(52)=1 - use the user's multi-channel optimisation and VEGAS grid 

files. 

ACSET(52)=2 - use the default multi-channel optimisation and user's VEGAS 
grid files. 

ACSET(52)=3 - use the default multi-channel optimisation and VEGAS grid 
files; read the user maximal weight file. 

• ACSET(53): Maximum weight search 

Mode for the maximum weight search needed for unweighting procedure: 
ACSET(53)=0 - no, use the provided files containing maximal weights. 
ACSET(53)=1 - use the provided files for max. weights and re-calculate the 
max. weights using the stored 100 events with the highest weight. 
ACSET(53)=2 - perform the search and give the new wtmax_xx_new.dat file; 
the switch is equivalent to generation of weighted events. 

• ACSET(54): Maximum weight choice 

Use the a-cutoff maximal weight wt^ ax or the overall maximal weight wt max 
found in training (see Section 4.2.3 for the explanation on these two options). 
ACSET(54)=0 - use the wt^ ax weight. 
ACSET(54)=1 - use the wt max weight. 



6.3 How to prepare data-files for the non-default setup 

The following actions are possible, to recover better efficiency of the generator 
modules with the non-default settings: 

• The user wants to generate events using different parton density function 
sets and/or different coupling values (e.g. AcerMC third order a s (Q 2 ) in- 
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stead of the first order one): 

It should suffice to set the the switch ACSET(53)=1, which signals Ac- 
erMC to re-calculate the wt max and wt^ ax using the 100 events stored in 
the file 

acermc_dat/grids/vtmaxA_xxYYY.dat. The coupling and parton density 
functions values should not change significantly the process topology but 
affect foremost the overall scale of the event weights; thus, the stored hun- 
dred events should still remain the ones with the highest weights and the 
re-calculated approximate estimates of the highest weight should be accu- 
rate enough. 

In case the user is not confident in the obtained result, the new maxi- 
mal weight estimation can be initiated by setting the switch ACSET(53)=2, 
which will result in generation of weighted events. The number of gener- 
ated events is determined by the usual NEVENT in run . card. At the end 
of the run AcerMC will produce a file called wtmax_xx_new.dat, with xx 
specifying the process number. The user should then start the generation of 
unweighted events with the setting ACSET(52)=3 and linking(renaming) the 
new file to wtmax_xx_usr.dat, with xx denoting the process number (e.g. 
wtmax_01_usr . dat). 

• The user wants to generate events using different values of particle/boson 
masses or other significant changes of the parameters apart from the centre- 
of-mass energy and/or m z o/ 7 * cutoff value for processes 5-8: 

In this case the user should re-train the VEGAS grid since the process 
topology is assumed to undergo minor changes. This is done by setting the 
switch ACSET(50)=2 or ACSET(50)=3; in the first case AcerMC starts with 
an untrained grid and in the second one it starts modifying the existing grid 
provided for the process at the selected hard process scale. In general the 
second option should be preferable since the topology should still be close 
to the pre-trained one. AcerMC again produces weighted events and at the 
end of the run outputs a file grid_xx_new. veg. The number of generated 
events is determined by the usual NEVENT in run. card. As in the previous 
case, the user should re-name the file to grid_xx_usr .veg, re-set the switch 
to ACSET(50)=0 and repeat the maximal weight search procedure described 
above, by setting the switch ACSET(53)=2 etc.. When the maximum weight 
search is completed the user switch ACSET(52)=2, which will cause Ac- 
erMC to read the wtmax_xx_usr.dat as well as grid_xx_usr .veg files and 
produce unweighted events with the new setup. 

• The user wants to generate events at a different mzo/ 7 * cutoff value and/or 
different centre- of-mass energy \fs: 

When the user changes at least one of these two parameters the event 
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topology is significantly changed as well as the contributions from different 
kinematic channels. The user should thus start with a new multi-channel 
optimisation by setting the mode switch ACSET(50)=1 and start an Ac- 
erMC run. The number of generated events is determined by the usual 
NEVENT in run. card. At the end of the run AcerMC will produce a file 
chanwt_xx_new. dat which should be renamed/linked to chanwt_xx_usr . dat. 
The user should then set the switch ACSET(52)=1 and first put ACSET(50)=2 
and and repeat the VEGAS grid training as described above and consequently 
ACSET(50)=0 and ACSET(53)=2 to perform the maximum weight search. Af- 
ter obtaining all three user files the ACSET(53) should again be put back to 
ACSET(53)=0 and a normal run should be started; the switch ACSET(52)=1 
will in this case force AcerMC to read all three user files and produce 
unweighted events. 

At the first look procedure for listed action scenaria might seem a bit complex 
but should after a few trials and errors become a straightforward routine; it 
is expected that the vast majority of users would have to deal with at most 
the first scenario. 
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6.4 Details on the interface to PYTHIA 6.2 

The AcerMC interface to Pythia 6 . 2 is implemented close to the new stan- 
dard specified at the Les Houches workshop 2001 [3]. The full description 
of the standard can be found in the PYTHIA 6.2 manual ([1]). In addition to 
the UPINIT and UPEVNT routines the file acermc_src/ interf ace/pythia_ac . f 
provides links between a list of AcerMC routines and the corresponding 
PYTHIA ones, as e.g. the (pseudo-) random number generator, a s and ojqed cal- 
culations as well as a series of routines that re- write the AcerMC event output 
to the required PYTHIA format. Using this strategy, the native AcerMC code 
is completely de-coupled from the linked hadronisation library (at the moment 
PYTHIA/HERWIG) and new interfaces can thus easily be added. The special Ac- 
erMC requirement is the call to the ACFINAL subroutine at the end of the run 
which signals the AcerMC to close the various I/O files and produce the final 
output. An example of the implementation of the PYTHIA/AcerMC interface 
can be found in the provided demo_py . f . The PYTHIA code is unmodified apart 
from making a small modification in PYIMIT routine: 

CALL UPINIT (1) 
..parameter inita lisation . . 
CALL UPINITC2) 
..process initialisation.. 

since the user-supplied processes in this new interface are not allowed to (re- 
)estimate maximal weights (as e.g. the native PYTHIA processes do). In the 
original code the call to UPINIT is set before the PYTHIA parameters and 
functions (e.g. PYALPS for a s (Q 2 ) calculation) are initialised with the user 
settings^]. 

The user can add the most recent PYTHIA library without other modifications 
but for the two lines of code in PYINIT routine as described above (the dummy 
routines UPINIT, UPEVNT, STRUCTM , STRUCTP and PDFSET however have to 
be removed from the code for an external process to work and to activate the 
PDFLIB interface). 

By setting ACER=1 user decides to generate hard process from AcerMC li- 
brary. Modeling of ISR/FSR shower, hadronisation and decays are generated 
by PYTHIA generator. All steering parameters, relevant for these steps of full 
event generation remain the same as in standard PYTHIA execution. 

By setting ACER=0 user decides to generate standard PYTHIA process. The 
simple example how to generate ttH production process within AcerMC 
framework is provided in demo_py.f. 



This was however possible in the old PYTHIA interface 
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6.5 Details on the interface to HERWIG 6.3 

Interfacing the AcerMC to HERWIG 6.3 has proved to be more of a prob- 
lemf Tr |, since the interfaced version does not comply with the Les Houches 
standard [3]. The interface routines are nevertheless written in the way to 
mimic the Les Houches description; in the HERWIG 6 . 3 interface the UPINIT 
routine however has to be called in two steps in order to enable the user to 
add the changes of the HERWIG 6 . 3 internal settings: 

CALL UPINIT (1) 
..user values.. 
CALL UPINIT (2) 

The UPEVNT subroutine in the HERWIG case replaces the native HWEPRO sub- 
routine as an iterface to AcerMC processes; in addition, it enables a user 
to activate more than one process in the same run, which is absent from the 
native HWEPRO implementation. As in the PYTHIA implementation all the inter- 
face subroutines needed for communication between HERWIG and AcerMC are 
stored in acermc_src/interf ace/herwig_ac . f . Some minor modifications of 
the original HERWIG 6 . 3 code were albeit necessary: 

• IMPLICIT NONE was commented out in the herwig6301 . inc file; this was 
needed since the AcerMC code is written with the implicit IMPLICIT 
DOUBLE PRECISION(A-H,0-Z). 

• The HERWIG interface to PDFLIB was changed in order to link it to the latest 
PFFLIB804 implementation which has a different PDFSET syntax. 

In principle the implemented changes should be very easy and transparent 
for the transfer into new HERWIG releases; it is also anticipated that a future 
HERWIG version will already have the interface to PDFLIB updated accordingly, 
and interface to Les Houches standard implemented. An example of the use 
of AcerMC/HERWIG interface is provided in the file demoJiw.f. 

By setting ACER=1 the user decides to generate a hard process from AcerMC 
library. Modeling of ISR/FSR shower, hadronisation and decays are generated 
by HERWIG generator. All steering parameters, relevant for these steps of full 
event generation remain valid as for the standard HERWIG execution. 

By setting ACER=0 the user decides to generate standard HERWIG processes. The 
simple example how to generate Zbb production process within the AcerMC 



7 This interface is temporary and we plan to move to HERWIG 6.4 version in the 
next future 
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framework is provided in demoJiw. f. 



The output logs of the run are produced in the directory prod, the acermc . out 
file containing the AcerMC specific information and the outputs pythia. out 
and/or herwig.out listing the outputs of the respective supervising genera- 
tors. The information about the input values of the steering files is stored in 
run. out in order to facilitate the event generation 'bookkeeping'. The sample 
outputs are given in Appendix C. 

6.6 Definition of the energy scale 

A few different values of scale Q 2 used in the evolution of parton density 
functions as well as the running couplings a s (Q 2 ) and q;qed(Q 2 ) can be set 
by the switch ACSET(2) (remember that the factorisation and renormalisation 
scales are assumed to be equal in AcerMC). Note that the correct value of 
the scale to be used for certain processes is in principle not known; what was 
implemented in AcerMC are the most probable/usual choices on the market; 
in measurements the best value will have to be determined by data analysis. 

• Processes 1, 2, 9: 
ACSET(2): (D=l) 

1 - Q 2 = s 

2 - Q 2 = E (Ph 2 + = < m 2 T > 
3-Q 2 = E(^ 2 )/4 =<p 2 T > 

4 - Q 2 = (m t + m H /2) 2 , m H = 120 GeV/c 2 

• Processes 3 — > 4: 
ACSET(2): (D=l) 

1 - Q 2 = M 2 W 

2 - Q 2 = s*^, where q = b,t 
Z-Q 2 = M 2 V + P T 2 V 

4 - Q 2 = 0.5 • (s* w + s* q9 ) + (pf) 2 , where q = b,t 

• Processes 5^8: 
ACSET(2): (D=l) 

1 - Q 2 = M§ is the only setting implemented at the moment. 

6. 7 Installation procedure 

The installation requires availability of the CERNLIB fortran library. 

• Ungzip and untar distribution file. 

• Go to the directory . /pythia_lib/src_pythia; type make install. It will 
compile the sources and install the created library libpythia. a to directory 
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pythia_lib. 

Go to the directory . /herwig_lib/src_herwig; type make install. It will 
compile the sources and install the created library libherwig. a to directory 
herwig_lib. 

Go to the directory . /helas_lib/src_helas; type make install. It will 
compile the sources and install the created library libhelas . a to directory 
helas_lib. 

Go to the directory . /acermc_src; type make install. It will the com- 
pile sources and install the necessary object files and the created library 
libacermc.a to the directory . /acermc_lib. 

In the main directory type make demo_py or make demo_hw. It will com- 
pile demo_py . f or demoJiw . f and produce the executables demo_py . exe or 
demoJiw . exe depending on the selected option. 

To execute the programs type run demo_py or run demoJiw. The scritps 
will change directory to prod and create respective links to data directories 
there. The execution will also be performed there. All input files should be 
accessible/routed from directory prod, the output files will also be produced 
in that directory. 
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7 Outlook and conclusions 

In this paper we presented the AcerMC Monte Carlo Event Generator, based 
on the library of the matrix-element-based generators and interfaces to the 
universal event generators PYTHIA 6 . 2 and HERWIG 6 . 3. The interfaces are 
based on the standard proposed in [3]. 

The presented library fulfills the following goals: 

• It gives a possibility to generate the few Standard Model background pro- 
cesses which were recognised as very dangerous for the searches for the 
New Physics at LHC, and generation of which was either unavailable or not 
straightforward so far. 

• Although the hard process event is generated with matrix-element-based 
generator, the provided interface allows to complete event generation with 
initial and final state radiation, multiple interaction, hadronisation, frag- 
mentation and decays, using implementation of either PYTHIA 6 . 2 or HERWIG 
6.3. 

• These interfaces can be also used for studying systematic differences between 
PYTHIA 6 . 2 or HERWIG 6 . 3 predictions for the underlying QCD processes. 

The complete list of the native AcerMC processes implemented so far is: 
99, QQ ttbb, qq -> £u)bb, gg,qq -> Z/Y{~^ M)Vb, qq -> W(-> lv)tt, 

99, QQ —> ££,w,bb)tt and the complete EW gg — > (Z/W/j* —>)ttbb. 

We plan to extend this not too exhaustive, but very much demanded list of 
processes, in the near future. 

Several improvements of the existing Monte Carlo algorithms/programs have 
been developed in the process of this work. Let us make short list of the most 
interesting ones: (l)The colour flow information has been obtained after some 
modification of MADGRAPH package; (2) The power of the multi-channel opti- 
misation was enhanced by using the modified ac-VEGAS package. We believe 
that the modification in the VEGAS code represents a very powerful extension of 
this package; (3) The additions and extensions to the available (multi-channel) 
phase space algorithms (e.g. Breit-Wigner function with s-dependent width) 
lead to substantial improvement of the unweighting efficiency; Figs. 3, 5 and 6 
illustrate the improvements achieved in the generation efficiency. 

Having all these different production processes implemented in the consistent 
framework, which can be also directly used for generating standard processes 
available in either PYTHIA 6 . 2 or HERWIG 6 . 3 Monte Carlo, represents very 
convenient environment for several phenomenological studies dedicated to the 
LHC physics. Such frame was not available to our knowledge so far. We hope 
that it can serve as an interesting example or even a framework. This way 
some tools for discussing the ambiguities due to QCD effects are collected, 
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however the necessary discussion for the appropriate uncertainties is still not 
exhausted. Nevertheless some discussions using this tool can be already found 
in [7], [8], [9], [12]. 
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A Feynman Diagrams 



The 38+7 Feynman diagrams contributing to the gg, qq — > ttbb production. 
Only four flavours are included for incoming quarks. Contribution of the in- 
coming b-quarks could be excluded from the calculations thanks to very high 
suppression induced by either the parton density functions and/or CKM ma- 
trix elements. 



Diagrams by MadGraph Diagrams by MadGraph 




Diagrams by MadGraph Diagrams by MadGraph 




graph 37 graph 38 graph 7 



Fig. A.l. The Feynman diagrams for the processes gg,qq — > ttbb. 
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The Feynman diagrams contributing to the qq — * W(— > £z/)66 and gg — > VF(— ■> 
£z/)tt matrix element are just two t-channel diagrams with fermion exchange 
and double conversion into an off-shell W boson and a virtual gluon; the W 
boson subsequently decays leptonically into Iv and the gluon splits into a bb 
pair or ti pair respectively, (c.f. Figure A. 2). 



Diagrams by MadGraph 




graph 1 graph 2 



Fig. A. 2. The Feynman diagrams for the process qq — > Wbb — > e + ^ e 66. T/ie same 
se£ zs used /or gg — > — ► e + z^ e it process, with b-quarks replaced by top-quarks. 
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The Feynman diagrams contributing to the gg, qq — > Z/ r y*{f f , vv)bb produc- 
tion are shown in Figure A. 3. The dominant contribution comes from the (2) 
and (6) configurations for the processes with gg initial state and the double 
conversion configuration (2), (4) for the ones with qq initial state. The same 
set of Feynman diagrams is used for the gg, qq — > Z/j*(ff, uu)tt process, with 
b-quarks being replaced by the top-quarks. If the Zj^f{-^ bb) decay mode is 
simulated, it represents only subset of the EW production of ttbb final state. 
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The complete set of th the Feynman diagrams contributing to the full electro- 
weak gg — > (Z/W/Y —>)bbti production mediated by exchange of the Z/W/j* 
bosons is shown in Fig. A. 4. 



Diagrams by MadGraph 



Diagrams by MadGraph 
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Diagrams by MadGtaph 




A. 4. The Feynman diagrams for the processes gg — > (Z/W/j* 



63 



B Example input files 



B.l File run. card 



c 

C STEERING FILE FOR ACERMC (l.OBETA) - BASIC SETTINGS 

c 

C==== TURN ON FFKEY STEERING FILE (DEBUG) 
LIST 

C CMS/ACSET(1) 

C Specify the centre-of-mass energy in GeV 
CMS 14000.0 

C==== ACER 
C Use AcerMC code 
C ACER=1 - yes 
C ACER=0 - no 

ACER 1 

C==== PRDCESS/IACSET ARRAY 

C Specify the process to generate. The available AcerMC processes are: 
C 

C 1) g + g -> t t' b b~ (MG) 

C 2) q + q~ -> t t~ b b~ (MG) 

C 3) q + q~ -> (W-» 1 nu_l b b~ (MG) 

C 4) q + q~ -> (W-» 1 nu_l t t" (MG) 

C 5) g + g -> (Z0->) 1 1" b b" (MG) 

C 6) q + q~ -> (Z0->) 1 1" b b" (MG) 

C 7) g + g -> (Z0->) f f" t t" (MG) 

C 8) q + q~ -> (Z0->) f f" t t" (MG) 

C 9) g + g -> (ZO/W/gamma->)t t~ b b~ (MG) 

C 

C In case ACER=0 the native Pythia/Herwig conventions should be used 
PROCESS 9 

C==== HAD 

C Control of hadronization/f ragmentation/ISR/FSR switches: 

C HAD=0 - switch off radiation in initial and final state, multinteraction and 
C hadronization 

C HAD=1 - switch off radiation in final state and hadronization 
C HAD=2 - switch off hadronization 
C HAD=3 - full treatment 

C HAD=4 - switch off radiation in initial state, multiinteraction and hadronization 
C HAD=5 - switch off radiation in initial state and multiinteraction 

HAD 3 

C==== PDF-GR0UP/ACSET(6) 

C Choose a PDF group according to PDFLIB804 naming scheme 
PDFGROUP 4 

C==== PDF-SET/ACSET(7) 
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C Choose a PDF set according to PDFLIB804 naming scheme 

C PDFSET=55 in not (yet) in the manual but represents CTEQ5L parametrised set 
PDFSET 55 

C==== RSEED 

C Choose the random seed for random generator initialisation 
RSEED 611123 

C==== NEVENT/ ACSET (51) 

C Specify the number of events to generate 
NEVENT 100 

c 

END 



B.2 File acermc.card 



c 

C STEERING FILE FOR ACERMC (1.0BETA) - ACERMC SETTINGS 

c 

C==== TURN ON FFKEY STEERING FILE (DEBUG) 
LIST 

C THE AcerMC EVENT SETTINGS 

C==== SCALE/ACSET(2) 

C Choose the Q~2 scale for the active AcerMC process. 

C The implemented values differ for various processes, please look into the manual 
C for details 

ACSET2 4 

C FERMI ON/ ACSET (3) 

C The flavour of the final state leptons produced in W or Z decays of AcerMC 
C processes 3-6. The Pythia/PDG naming convention is used: 
C FERMI0N=11 - electron 
C FERMI0N=13 - muon 
C FERMI0N=15 - tau 

C FERMI0N=12, 14, 16 - neutrinos, a single flavour is generated and the cross-section is 
C calculated accordingly. For a neutrino final state the user should 

C generate all three flavours! 

C FERMI 0N=5 - b-quark final state 

C The last two settings ACSET (3) =12, 14, 16 ,5 work only for processes 7 and 8! 
ACSET3 5 

C==== Z/GAMMA/ACSET (4) 

C Use the full Z/gamma* propagator in AcerMC processes 5-6. 

C ZGAMMA=0 - only Z propagator 

C ZGAMMA=1 - full Z/gamma* propagator 

ACSET4 1 

C Z/G CUT/ACSET(5) 
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C Cutoff value on the invariant mass m_Z/gamma* in GeV when ZGAMMA=1 . 
C Note that the provided data files exist only for values of 
C ZGCUT=2,5, 10, 15,30,60, 120,270,300,500 GeV which should satisfy most 
C users. In case a different value is set the user has also to provide 
C the user data files for the run. 

ACSET5 10.0 

C THE AcerMC ADVANCED SWITCHES 

C ALPHA_S/ACSET(8) 

C Use the alpha_s provided by the linked generator (Pythia/Herwig) or the 
C one provided by AcerMC 

C ALPHAS=0 - use the linked generator's alpha_s 
C ALPHAS=1 - use the AcerMC s alpha_s (one loop calculation) 
C ALPHAS=2 - use the AcerMC s alpha_s (three loop calculation) 
ACSET8 

C==== LAMBDA_S/ACSET(9) 

C Specify the value of lambda" (nf =5) _MSbar for AcerMC alpha_s calculation 

C LAMBDAS=-1 - the value is taken from PDFLIB804 for the corresponding 

C parton density function set 

C LAMBDAS>0 - the provided value is taken 

ACSET9 -1. 

C==== ALPHA_EM/ACSET (10) 

C Use the alpha_QED provided by the linked generator (Pythia/Herwig) or the 
C one provided by AcerMC 

C ALPHAEM=0 - use the linked generator's alpha_QED 
C ALPHAEM=1 - use the AcerMC 's alpha_QED 
ACSET10 

C==== ALPHA_EM(0)/ACSET(11) 

C Specify the value of alpha_QED(0) for AcerMC alpha_QED calculation 
C ALPHAEM0=- 1 - the default AcerMC value is used 
C ALPHAEMO0 - the provided value is taken 

ACSET11 -1. 

C==== TOP S-L/ACSET(12) 

C Specify the decay mode of WW pair produced by top decays in AcerMC 
C processes 1,2,4,7,8 and 9: 

C TOPDEC=0 - both W bosons decay according to Pythia/Herwig switches 

C T0PDEC=1 - one W decays into electron + nu and the other one hadronically 

C T0PDEC=2 - one W decays into muon + nu and the other one hadronically 

C T0PDEC=3 - one W decays into tau + nu and the other one hadronically 

C T0PDEC=4 - one W decays into el or mu + nu and the other one hadronically 

C T0PDEC=5 - one W decays into el or mu + nu and the other one hadronically, the 

C W decaying leptonically has the same charge as the primary W; the decay 

C mode makes sense only for AcerMC processes 4! 

C When T0PDEO0 the output cross-section is ALREADY MULTIPLIED by the corresponding 
C branching ratio! (Courtesy of AcerMC authors) 

ACSET12 

C THE AcerMC TRAINING SETUP AND UNWEIGHTING TREATMENT 

C==== M0DE/ACSET(50) 

C Specify the AcerMC training mode: 

C MODE=0 - normal run, generate unweighted events 

C M0DE=1 - perform multi-channel optimisation. 

C M0DE=2 - perform VEGAS grid training. 
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C M0DE=3 - perform VEGAS grid training as M0DE=2 but does this by updating a provided grid 
ACSET50 

C==== USER/ACSET(52) 

C Use the data files provided by user 
C USER=0 - no, use internal files 

C USER=1 - use the user's multi-channel optimisation and VEGAS grid files 

C USER=2 - use the default multi-channel optimisation and user's VEGAS grid files 

C USER=3 - use the default multi-channel optimisation and VEGAS grid files, read the user 

C maximal weight file. 

ACSET52 

C==== MAXFIND/ACSET(53) 

C Search for the maximum weight needed for event unweighting 
C MAXFIND=0 - no, use the provided file for max. weights 

C MAXFIND=1 - use the provided file for max. weights, re-calculate the max. weights using 
C the stored 100 highest events 

C MAXFIND=2 - perform the search and give the wtmax file, equivalent to generation of 
C weighted events 

ACSET53 1 

C==== EPSIL0N/ACSET(54) 

C Use the epsilon maximal weight or the overall maximal weight found in training (see the 

C manual for the difference) 

C EPSIL0N=0 - use the epsilon max. weight 

C EPSIL0N=1 - use the overall maximal weight 

ACSET54 

c 

END 
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C Example output files 

C.l File acermc.out 



-+++=+= 

: I :-: 1 = 

-; 1 = 



-r/.x]] 
-/] i 
_i] 



Ml. 
1 = 1; 



• j%3]:, 

-lyx/xm 

_ I xx'/,xxxx'/,7,'/,'/,'/,+ ' 
: ~] '/,xxxx] xx'/,x_ , +_x_% ' 
-__l Ix| |xx]+]]]+]]]]x|xxx] 1 
-+r/„xxxx] ]]]+]]++]+] x] I >- .;..;.:_/' 

-+x]]]] l+]+]+]++]++]] + l + |] + |]]+- 
., | x] +] + I I =++]+]=++++++= I ++]=+] l~ - 

] ;x]] | | =++++++++ | ] ; I ++++++++++=+] =] ; 

]+] |]++++++++| | | |==| : ; : =|==++==+; I , : 

I +>+++] +]]| + lll = l = ;l;;:|; = : 1==;::.;;: 

-++]+++:+ |x+: : I |=| I :=:::::;;::::::::;: ; 
I x+ 1+1]++:, -..:::! = ; = :-.:.:-::. 



J J J ■>>>>>> • 



-I;:, 



|;:|:::. - . : . - : I ; I I : : . - 
|==|:;:. :-:-:. 

-::::. I I i === :::.:.::: 
-:=;;l+ll=;==l=::::: 

; ; ; : l-l+=+l + l==; : ' 
. 1=1 l=+l=:=, 

-+;,□! ;; = ' 
- :l; - 



40000L, 
~?##0L 
dO 40, 

J01 -#w 

_00yyyW0L 
##!!!! !#0; 



. aaaa aa 
_W0#V9N0#& 
#0' ?#W 

:0f 

-OA 



. aaaa ; 
dO#V*N#0, 
##~ -#0; 
:000###00001 
-OA 



aaaa, _aaa, 
0##0LW0@4#<3 ; 

j##9 

j#l 



lOOOOi j000& .a00000L#0 
-000A _0001- _d0! ' 
00j#; JO I 01 d0> 
00 4# I 01 1 01 00 

oo ?#oa' |oi #0 



-400 
40 



00 HH< 101 jOL 



adOLa, ajOAa 4#Aaa_aj#0' ?OLaa_aaaOL aaJOLaaa, _a00aa _ajOLa *0Aaa_aad 



HHHRH1 HHHHH 



1 9##009 ! 



'9NW009! ! 



HHHHRHHH1 : HHHHH 7HHHRH ?!##00P!' 



AcerMC 1.0 (February 2002), B. P. Kersevan, E. Richter-Was 



< ACTIVATED PROCESSES > 



1) 


g 


+ 


g 


-> 


t t~ b b" (MG) 




OFF 


2) 


q 


+ 


q" 


-> 


t t~ b b" (MG) 




OFF 


3) 


q 


+ 


q" 


-> 


(W->) 1 nu_l b b" 


(MG) 


OFF 


4) 


q 


+ 


q" 


-> 


(W->) 1 nu_l t t" 


(MG) 


OFF 


5) 


g 


+ 


g 


-> 


(Z0->) 1 1" b b~ 


(MG) 


OFF 


6) 


q 


+ 


q" 


-> 


(Z0->) 1 1" b b~ 


(MG) 


OFF 


7) 


g 


+ 


g 


-> 


(Z0->) f f" t t~ 


(MG) 


OFF 


8) 


q 


+ 


q" 


-> 


(Z0->) f f" t t~ 


(MG) 


OFF 


9) 


g 


+ 


g 


-> 


(Z0/W/gamma->) t t 


- b b~ 


(MG) ON 



■< ACERMC SETTINGS > 
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C.M.S ENERGY = 
SCALE CHOICE = 
ACERMC ALPHA_QCD 
LAMBDA (5) _MS 
ACERMC ALPHA_QED 
ALPHA_QED(0) 
TOP->W S-L DECAY 



14000.00 
4 


-1.000000 


-1.000000 




[ACSET(l)] 

[ACSET(2)] 

[ACSET(8)] 

[ACSETO)] 

[ACSET(IO)] 

[ACSET(ll)] 

[ACSET(12)] 



OPTIMIZATION 
OPTIM. STEPS 
USER FILES 
MAX. SEARCH 
EPSILON CUTOFF 




100 

1 




[ACSET(50)] 
[ACSET(51)] 
[ACSET(52)] 
[ACSET(53)] 
[ACSET(54)] 



-< APPROXIMATE MAXIMUM WEIGHT ESTIMATION >- 



NEW MAXIMUM WEIGHT (MB) = 0.991294E-08 
NEW EPSILON WEIGHT (MB) = 0.865417E-08 



-> APPROXIMATE MAXIMUM WEIGHT ESTIMATION <- 



■< FINALIZATION FOR PROCESS: 9 > 

< WEIGHT SURVEY > 

< TOTAL STATISTICS > 



CROSS-SECTION ESTIMATE = 
+/- 

VARIANCE ESTIMATE = 
+/- 



0.877415E+00 PB 
0.378414E-01 PB 
0.143197E-02 PB~2 
0.107049E-03 PB~2 



MAXIMUM WEIGHT = 0.656660E-08 
NO. WEIGHTS NE = 1009 
NO. WEIGHTS EQ = 
NO. WEIGHTS LT = 

EFFICIENCY FOR ALL WEIGHTS = 13.362 '/, 

EFFICIENCY FOR NONZERO WEIGHTS = 13.362 '/„ 

NO. WEIGHTS ABOVE EPSILON-CUT = 

> WEIGHT SURVEY < 
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C.2 File pythia.out 



** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 



Welcome to the Lund Monte Carlo ! 



lh 



hh 



11 



PPP 
P P 
PPP 
P 
P 



Y Y 
Y 
Y 
Y 



Y TTTTT H H III 
H H I 
HHHHH I 
H H I 
H H III 



A 
A A 
AAAAA 
A A 
A A 



This is PYTHIA version 6.203 
Last date of change: 13 Nov 2001 

Now is Jan 2000 at 0:00:00 

Disclaimer: this program comes 
without any guarantees. Beware 
of errors and use common sense 
when interpreting results. 

Copyright T. Sjostrand (2001) 



** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 

** ** 
** ** 

***************************************************************** 
1****************** PYINIT: initialization of PYTHIA routines ***************** 



** An archive of program versions and documentation is found on the web: 

** http : //www. thep. lu. se/~torbjorn/Pythia. html 

** 

** When you cite this program, currently the official reference is 

** T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Hrenna and 

** E. Norrbin, Computer Physics Commun. 135 (2001) 238. 

** The large manual is 

** T. Sjostrand, L. Lonnblad and S. Mrenna, LU TP 01-21 [hep-ph/0108264] . 

** Also remember that the program, to a large extent, represents original 

** physics research. Other publications of special relevance to your 

** studies may therefore deserve separate mention. 
** 

** Main author: Torbjorn Sjostrand; Department of Theoretical Physics 2, 

** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; 

** phone: + 46 - 46 - 222 48 16; e-mail: torbjornSthep. lu. se 

** SUSY author: Stephen Mrenna, Physics Department, UC Davis, 

** One Shields Avenue, Davis, CA 95616, USA; 

** phone: + 1 - 530 - 752 - 2661; e-mail: mrennaOphysics.ucdavis.edu 



PYTHIA will be initialized for p+ on p+ user configuration 
with 7000.000 GeV on 7000.000 GeV beam energies 

corresponding to 14000.000 GeV center-of-mass energy 



******** PYMAXI: summary of differential cross-section maximum search ******** 



I ISUB Subprocess name 
I 



I I 
I Maximum value I 
I I 
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User process 669 



8.6542E-09 I 
I 



********************** PYINIT: initialization completed ********************** 



Event listing (summary) 



I particle/jet KS 


KF 


orig 


P- 


.X 


P- 


-y 


P- 


_z 


E 


m 


1 !p+! 




21 


2212 








000 





000 


7000 


000 


7000 


000 





938 


2 !p+! 




21 


2212 








000 





000- 


-7000 


000 


7000 


000 





938 


3 !g! 




21 


21 


1 


1 


144 





067 


1425 


499 


1425 


499 





000 


4 !g! 




21 


21 


2 


1 


111 





247- 


-1179 


739 


1179 


739 





000 


5 !g! 




21 


21 


3 


1 


144 





067 


1425 


499 


1425 


499 





000 


6 !g! 




21 


21 


4 


1 


111 





247- 


-1179 


739 


1179 


739 





000 


7 !t! 




21 


6 





-24 


631 


-28 


056 


43 


153 


184 


068 


175 


000 


8 !tbar! 




21 


-6 





80 


035 


-23 


611 


1084 


153 


1101 


351 


175 


000 


9 !b! 




21 


5 





-30 


160 


24 


308 


217 


096 


220 


578 


4 


800 


10 !bbar! 




21 


-5 





-22 


988 


27 


673- 


-1098 


642 


1099 


241 


4 


800 


11 t 


A 


2 


6 


7 


-24 


631 


-28 


056 


43 


153 


184 


068 


175 


000 


12 uu_l 


V 


1 


2203 


2 


-0 


215 


-0 


195- 


-5351 


371 


5351 


371 
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13 tbar 


A 


2 


-6 


8 


80 


035 


-23 


611 


1084 


153 


1101 


351 


175 


000 


14 u 


V 


1 


2 


1 


-0 


494 


-0 


216 


1493 


138 


1493 


138 





330 


15 b 


A 


2 


5 


9 


-30 


160 


24 


308 


217 


096 


220 


578 


4 


800 


16 ud_0 


V 


1 


2101 


1 


-0 


651 





149 


4081 


363 


4081 


363 





579 


17 bbar 


A 


2 


-5 


10 


-22 


988 


27 


673- 


-1098 


642 


1099 


241 


4 


800 


18 d 


V 


1 


1 


2 


-0 


896 


-0 


052 


-468 


889 


468 


890 
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sum: 2.00 0.00 0.00 0.00 14000.00 14000.00 

1********* PYSTAT: Statistics on Number of Events and Cross-sections ********* 



I 
I 
I 

I 

I 

I N:o Type 
I 



Subprocess 



Number of points 



Generated 



I 
I 
I 

1 

I 

Tried I 
I 



Sigma 
(mb) 



I I 

I All included subprocesses I 

I 4 User process 669 I 

I I 



100 
100 



I I 

1009 I 8.774E-10 I 

1009 I 8.774E-10 I 

I I 



********* Fraction of events that fail fragmentation cuts = 0.00000 ********* 
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C.3 File herwig.out 



HERWIG 6.301 9 July 2001 

Please reference: G. Marcheslnl, B.R. Webber, 

G.Abbiendi, I.G.Knowles, M.H.Seymour & L.Stanco 
Computer Physics Communications 67 (1992) 465 
and 

G.Corcella, I.G.Knowles, G.Marchesini, S.Moretti, 

K.Odagiri, P.Richardson, M.H.Seymour & B.R.Webber, 
JHEP 0101 (2001) 010 

INPUT CONDITIONS FOR THIS RUN 

BEAM 1 (P ) MOM. = 7000.00 

BEAM 2 (P ) MOM . = 7000.00 

PROCESS CODE (IPROC) = 90660 

NUMBER OF FLAVOURS = 6 

STRUCTURE FUNCTION SET = 8 

AZIM SPIN CORRELATIONS = T 

AZIM SOFT CORRELATIONS = T 

QCD LAMBDA (GEV) = 0.1460 

DOWN QUARK MASS = 0.3200 

UP QUARK MASS = 0.3200 

STRANGE QUARK MASS = 0.5000 

CHARMED QUARK MASS = 1.5500 

BOTTOM QUARK MASS = 4.8000 

TOP QUARK MASS = 175.0000 

GLUON EFFECTIVE MASS = 0.7500 

EXTRA SHOWER CUTOFF (Q)= 0.4800 

EXTRA SHOWER CUTOFF (G)= 0.1000 

PHOTON SHOWER CUTOFF = 0.4000 

CLUSTER MASS PARAMETER = 3.3500 

SPACELIKE EVOLN CUTOFF = 2.5000 

INTRINSIC P-TRAN (RMS) = 0.0000 
MIN MTM FRAC FOR ISR =1.0000E-04 
1-MAX MTM FRAC FOR ISR =1.0000E-06 

NO EVENTS WILL BE WRITTEN TO DISK 

B_d: Delt-M/Gam =0.7000 Delt-Gam/2*Gam =0.0000 
B_s: Delt-M/Gam = 10.00 Delt-Gam/2*Gam =0.2000 

PDFLIB USED FOR BEAM 1: SET455 OF AcerMC 
PDFLIB USED FOR BEAM 2: SET455 OF AcerMC 



Checking consistency of particle properties 
Checking consistency of decay tables 



INPUT EVT WEIGHT = 6.9480E-03 
INPUT MAX WEIGHT = 6.9480E-03 

SUBROUTINE TIMEL CALLED BUT NOT LINKED. 
DUMMY TIMEL WILL BE USED. DELETE DUMMY 
AND LINK CERNLIB FOR CPU TIME REMAINING. 



EVENT 1: 7000.00 GEV/C P ON 7000.00 GEV/C P PROCESS: 90669 
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SEEDS : 



17673 & 



63565 STATUS : 



10 ERROR: 



WEIGHT: 8.2681E-04 



-INITIAL STATE- 



IHEP ID 


IDPDG 


1ST 


M01 


MO 2 


DAI 


DA2 


P-X 


P-Y 


P-Z 


ENERGY 


MASS 




v-x 


1 P 


2212 


101 

















00 





00 


7000 





7000 





0.94 





OOOE+OO 


2 P 


2212 


102 

















00 





00- 


•7000 





7000 





0.94 





OOOE+00 


3 CMF 





103 


1 


2 











00 





00 





014000 


014000.00 





OOOE+OO 


4 GLUON 


21 


111 


6 


10 





7 





00 





00 


378 


9 


378 


9 


0.00 





OOOE+00 


5 GLUON 


21 


112 


6 


8 





9 





00 





00 


-223 


6 


223 


6 


0.00 





OOOE+OO 


6 HARD 





110 


4 


5 


7 


10 





00 





00 


155 


3 


602 


5 


582.15 





OOOE+00 


7 TQRK 


6 


113 


6 


4 





8 


10 


38 


71 


83 


-18 


5 


190 


4 


175.00 





OOOE+OO 


8 TBAR 


-6 


114 


6 


7 





5 


37 


44 


-107 


58 


215 


3 


300 





175.00 





OOOE+00 


9 BQRK 


5 


114 


6 


5 





10 


-37 


78 


29 


50 


-67 


7 


83 


1 


4.80 





OOOE+OO 


10 BBAR 


-5 


114 


6 


9 





4 


-10 


05 


6 


25 


26 


2 


29 


1 


4.80 





OOOE+00 



EVENT 
SEEDS : 



2: 7000.00 GEV/C P 
1689786158 & 1462719137 



ON 7000.00 GEV/C P 
STATUS: 10 ERROR: 



PROCESS: 90669 
WEIGHT: 3.9180E-03 



-INITIAL STATE- 



IHEP 


ID 


IDPDG 


1ST 


M01 


MO 2 


DAI 


DA2 


P-X 


P-Y 


P-Z 


ENERGY 


MASS 




v-x 


1 


P 


2212 


101 

















00 





00 


7000 





7000.0 


0.94 





OOOE+OO 


2 


P 


2212 


102 

















00 





00- 


•7000 





7000.0 


0.94 





OOOE+00 


3 


CMF 





103 


1 


2 











00 





00 





014000.014000.00 





OOOE+OO 


4 


GLUON 


21 


111 


6 


8 





9 





00 





00 


506 





506.0 


0.00 





OOOE+00 


5 


GLUON 


21 


112 


6 


10 





7 





00 





00 


-142 


1 


142.1 


0.00 





OOOE+OO 


6 


HARD 





110 


4 


5 


7 


10 





00 





00 


363 


9 


648.1 


536 . 23 





OOOE+00 


7 


TQRK 


6 


113 


6 


5 





8 


-77 


38 


1 


14 


70 


1 


203.8 


175.00 





OOOE+OO 


8 


TBAR 


-6 


114 


6 


7 





4 


54 


35 


-16 


24 


326 


2 


374.5 


175.00 





OOOE+00 


9 


BQRK 


5 


114 


6 


4 





10 


-3 


39 


-5 


11 


11 


9 


14.2 


4.80 





OOOE+OO 


10 


BBAR 


-5 


114 


6 


9 





5 


26 


42 


20 


21 


-44 


3 


55.6 


4.80 





OOOE+00 



OUTPUT ON ELEMENTARY PROCESS 



N.B. NEGATIVE WEIGHTS NOT ALLOWED 



NUMBER OF EVENTS 
NUMBER OF WEIGHTS 
MEAN VALUE OF WGT 
RMS SPREAD IN WGT 
ACTUAL MAX WEIGHT 
ASSUMED MAX WEIGHT 



100 
1074 
6.9823E-04 
9.3117E-04 
5.1161E-03 
6.9480E-03 



PROCESS CODE IPROC 
CROSS SECTION (PB) 
ERROR IN C-S (PB) 
EFFICIENCY PERCENT 



90669 
0.6982 
2.8414E-02 
10.05 
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